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Instantons generate strong non-perturbative interactions between quarks. In the vacuum, these 
interactions lead to chiral symmetry breaking and generate constituent quark masses on the order 
of 300-400 MeV. The observation that the same forces also provide attraction in the scalar diquark 
channel leads to the prediction that cold quark matter is a color superconductor, with gaps as large 
as ~ 100 MeV. We provide a systematic treatment of color superconductivity in the instanton model. 
We show that the structure of the superconductor depends on the number of flavors. In the case 
0^ ' of two flavors, we verify the standard scenario and provide an improved calculation of the mass 

On , gap. For three flavors, we show that the ground state is color-flavor locked and calculate the chiral 

condensate in the high density phase. We show that as a function of the strange quark mass, there 
is a sharp transition between the two phases. Finally, we go beyond the mean-fleld approximation 
bJQ, and investigate the role of instanton-antiinstanton molecules, which - besides superconducting gap 

^ ' formation - provide a competitive mechanism for chiral restoration at finite density. 
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I. INTRODUCTION 



I> [ Studying Quantum Chromodynamics (QCD) at finite baryon density does not require a special motivation: this is, 
CO . after all, the traditional subject of nuclear physics. Nevertheless, the investigation of cold quark matter lay dormant 
' for some time and was only revived recently when it was realized that a number of exciting new phenomena can be 
^2 , predicted with some certainty. The first is that we expect the high density phase of QCD to be a color superconductor, 
with sizeable gaps on the order of 100 MeV around the phase transition ||l|,^. The structure of this phase depends 
, sensitively on the number of active flavors. For three or more flavors, color and flavor quantum numbers are locked 
^\ • and chiral symmetry is broken even at large chemical potential fj, In addition to that, it was argued that the 

second order chiral phase transition at finite temperature T and zero density is likely to turn first order at some 
critical density. This entails the existence of a tricritical point in the /i — T phase diagram, which would persist even 
if the light quarks are not massless 1^-01 . 
O-i. The idea that asymptotic freedom and the presence of a sharp Fermi surface imply that high density QCD should be 
^ ' a color superconductor goes back to the work of Frautschi, Barrois, Bailin and Love |P-p^. Color superconductivity 
has many features of the standard model, such as dynamical gauge symmetry breaking and the Higgs phenomenon. It 
is different from electroweak symmetry breaking in the standard model in the sense that the Higgs is composite. And 
it is different from models with compositeness (such as technicolor) in that it does not require strong interactions. 
• Color superconductivity takes place even in weak coupling. This, of course, is a consequence of the BCS instability. 
' Detailed numerical calculations of color superconducting gaps were carried out by Bailin and Love, who concluded 
that one-gluon exchange (OGE) induces gaps on the order of 1 MeV at several times nuclear matter density (at 
asymptotically large chemical potentials, however, magnetic gluon exchanges generate increasingly large gaps pd|). 
The main new feature pointed out in is that instanton-induced interactions can lead to substantially larger gaps, 
on the order of 100 MeV. Furthermore, it was realized that the phase structure of QCD at finite baryon density is 
very rich. Besides the dominant order parameter for the superconducting phase transition, which is a scalar-isoscalar 
color antitriplet diquark operator, many other forms are possible. 

Previous work mostly concentrated on two or three massless flavors and was based on the mean-field approximation 
(MFA). In the present work, we go beyond these approximations in several important respects. In the case of two flavors 
we replace schematic zero range interactions with the full momentum dependent instanton-induced interaction [p^|-p^ . 
We study the three-flavor case and show that the ground state exhibits color-flavor locking. We show that chiral 
symmetry is broken, calculate the chiral condensate, and also assess the effects of a finite strange quark mass. In 
order to go beyond the mean field approximation we study the role of instanton-antiinstanton clusters. In particular, 
we consider the competition between random instantons and clusters employing a statistical mechanics treatment of 
the partition function for the instanton liquid. Finally, we consider more speculative possibilities such as phases with 
diquark Bose condensation, and give a general discussion of the phase structure of QCD with different quark masses. 

Throughout the article we assume that instanton-induced effects are the predominant source of strong non- 
perturbative interactions in cold quark matter at small and moderate densities. This assumption is based on the 
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success of the instanton model at zero temperature and zero chemical potential, as well as at = and T 7^ 0, see 
for a review. The other major type of interaction that has been widely used is the (perturbative) OGE. While 
the latter should be prevailing at high densities, it encounters conceptual difficulties at low and moderate densities 
since the involved momentum transfers at the Fermi surface < and the running coupling constant of QCD 
might not be sufficiently small as required for a perturbative treatment. On the other hand, the Debye-screening of 
electric fields suppresses instanton effects at large densities; however, this suppression seems not to be effective below 
the chiral phase transition, as has been explicitely demonstrated in finitc-T lattice studies [T^Jl7| ] 

Further evidence for the importance of instantons in this context is related to the existence of fermion zero modes in 
the spectrum of the Dirac operator, which arise as a consequence of the axial anomaly: In a topologically non-trivial 
background field with topological charge +1 there is - for each flavor - a left-handed state that emerges out of the 
Dirac sea, and a right-handed one that moves from positive to negative energy. As a result, the axial charge is violated 
by 2Nf units. For Nf = 1 this immediately implies chiral condensation. For more than one flavor, chiral condensation 
is a collective effect. The quark condensate is determined by the number density of (almost) zero modes of the Dirac 
operator. These anomalously small eigenmodes can originate from the interaction of exact zero modes associated 
with isolated instantons and antiinstantons. The wave function of the condensate is the collective state built from 
instanton and antiinstanton zero modes. To check this mechanism for chiral symmetry breaking on the lattice has 
lately attracted appreciable attention: the results indeed support the suggested picture ||l^,0 . 

At finite baryon chemical potential the axial anomaly is connected with fermion zero modes in exactly the same 
way as in vacuum. The only difference is that the zero modes now correspond to extra states appearing at the Fermi 
surface, rather than the surface of the Dirac sea. As in the vacuum the effect of the zero modes can be represented 
as an effective (2A^j:)-quark interaction that operates near the Fermi surface. For two flavors, this interaction directly 
leads to the BCS instability, the formation of Cooper pairs and the appearance of a gap. For three and more flavors 
the instanton vertex does not directly support a Cooper pair; some of the chiral condensates have to be non-zero to 
close off external quark legs, reducing the A'^y^-body instanton interaction to a two-body one. Instanton-antiinstanton 
molecules, on the other hand, lead to an effective four quark operator for any number of flavors, which, if attractive, 
will trigger the formation of a gap. Nevertheless, direct instantons play an important role even for Nf > 3. In 
particular, instantons provide a novel mechanism for chiral condensation: A diquark-driven qq (chiral) condensate. 

The investigation of chiral symmetry restoration and color superconductivity at finite density should also be placed 
in a broader context, 6.17., including finite temperatures. Moreover, we would like to understand the phase structure 
as a function of parameters that we cannot control in the real world, such as the number of flavors and their masses. 
After all, the underlying mechanisms for the various transitions and the role of non-perturbative effects (such as 
instantons) in the different phases have to be clarified. 

At high temperature we expect to find a quark-gluon plasma phase in which chiral symmetry is restored, i.e., the 
density of (almost) zero modes has to vanish. This can be realized if the instanton liquid changes from a random 
ensemble of instantons and antiinstantons to a correlated system with finite clusters, e.g., instanton-antiinstanton 
(I-A) molecules. The formation of molecules and other correlated clusters was observed in numerical simulations 
of the instanton liquid where a number of consequences of this scenario were explored. On the lattice the 
disappearance of the quasi-zero modes in the vicinity of Tc is well established, and the formation of clusters has been 
observed Nevertheless, many details of the transition remain to be understood. In the case of many flavors 

the instanton calculations suggest a chirally restored vacuum state already for a fairly small number of flavors, 
around 5. Again, the transition is associated with the formation of correlated clusters. 

In this article we would like to understand the interplay of the three major phases that have been considered: (i) the 
hadronic (H) phase, with (strongly) broken chiral symmetry (ii) the color superconductor (CSC) phase, with broken 
color symmetry, and (iii) the quark-gluon plasma (QGP) phase. All three phases are associated with three specific 
instanton-induced interactions. Chiral symmetry breaking is caused by the strong qq attraction. The binding energy 
of the lightest baryon, the nucleon, is mostly associated with the qq interaction. The same interaction is responsible 
for superconductivity at large baryon density. Finally, quark exchanges between instantons and antiinstantons drive 
their pairing, which is expected to become the predominant feature in the instanton liquid as temperature increases 
(at any chemical potential). 

The structure of our paper is as follows. The first part, comprised of sects. II- VII contains a mean-field analysis 
of color superconductivity in finite-density QCD with 2 and 3 flavors. We begin with a brief introduction to the 
structure of the effective instanton-induced interaction in sect. In sect. [II, we study the physical effects of this 



interaction in different diquark channels at = 0. In sect. we discuss the modifications of the instanton-induced 
interactions at non-zero chemical potential. The interplay of chiral symmetry breaking and quark superconductivity 
in two-fiavor QCD is studied in sect. ^ This section employs the mean- field approximation but is based on the 
exact form of the interaction. Using a simplified version of the form factors we then consider different (qq) and (qq) 
condensates of increasing complexity: the three-flavor probl em i n the chiral limit (sect. VI) and the effect of flavor 
symmetry breaking due to a finite strange quark mass (sect. VII). 
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The second part of th e artic le (sect. VID and IX) addresses effects due to clustering, which are beyond the mean- field 
approximation. In sect. VIII we quantitatively discuss only one type of cluster, the instanton-antiinstanton molecule, 
which we believe is the most important cluster in the chirally restored phase. In sect. [X, we also discuss the role of 
correlations between quarks - in particular non-condensed diquarks and (the most obvious cluster of all!) nucleons 
in nuclear matter -, and comment onpossible experimental consequences for heavy-ion reactions and neutron stars. 
We summarize and conclude in sect. K. 



II. EFFECTIVE INSTANTON-INDUCED INTERACTIONS IN VACUUM 



A. Single-Instanton Interactions 

Our starting point is the euclidean QCD partition function 

Z = J Vi^Vi^^'VAexpi^SQCD) = J VAdet{p)exp{~S gauge) ■ (1) 

The main assumption of the instanton model is that the gauge field is saturated by classical (anti-)instanton solutions. 
If the instanton ensemble is sufficiently dilute, the gauge field can be approximated by a sum of individual instanton 
gauge potentials 

A = ^ ylfc . (2) 
keij 

Collective effects related to chiral symmetry breaking are generated through the low-momentum part of the fermion 
determinant. In particular, we will concentrate on the fermion determinant in a basis spanned by the zero modes of 
the individual instantons. Matrix elements of the Dirac operator in this basis are given by the overlap integrals 

Tja{z,u)= / d'^x (j)\{x - zj) ]P4>a{x-za)- (3) 



Here, zj and za denote the positions of the instanton and antiinstanton, and (pi^A the corresponding zero mode wave 
functions, which are solutions of the Dirac equation 

]Pi,AhAx)=^ . (4) 

where the covariant derivative pi, a includes the gauge potential of the (anti-) instanton / [A). Using the Dirac 
equation (^) and the sum ansatz (|^), Eq. (^) can be simplified by replacing the covariant derivative by an ordinary 
one. The overlap matrix element can also be viewed as the quark "hopping" amplitude from an instanton to an 
antiinstanton. 

To extract effective 2A'^/-quark interaction vertices, we follow the approach of Diakonov and Petrov who 
suggested to reintroduce free fermion fields according to 



1=1 1=1 



where in the two-flavor case 
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9+= J dQip^ [ Jd^xi>}{x)i ^hix-zi) J d'^y^iy - zi)z ^^/(y)] , 



(6) 



and the integrals are over the collective coordinates f2/ = {zj, pi,ui} (position, size and color orientation) of the 
instantons. The original zero mode determinant can be recovered by calculating a Green's function with Nf{N^+N_) 
external legs. In order to perform the integration over the centers of the instantons it is convenient to proceed to 
momentum space. This automatically induces a four-momentum conserving (5-functions at each vertex. 

An effective interaction is most easily derived by exponentiating the fermion terms. This is accomplished by applying 
an inverse Laplace transformation which gives the following partition function: 
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Z = const J di^di^Up+dp^ exp |-(A^+ + 1) log - (iV_ + 1) log 



13- 



c„ 



+ / d^x (V^i + f3+e+ + (3^6^) \ . (7) 



The integrations over (3± can be performed by the saddle point method, which becomes exact in the thermodynamic 
limit as the coefficients in the exponent are extensive quantities {N± = n±V4). For an equal number of instantons and 
antiinstantons, one may consider g = /3+ = /3_ as an effective fermion coupling. (3± are then eliminated through the 
final minimization of the free energy, leaving the total instanton density N/V = n+ + n_ as the physical parameter. 

In the remainder of this section we restrict ourselves to two fiavors. In this case, four quarks participate at each 
vertex, and the pertinent vertex operator 6± takes the form 

2 

Og^ = n dnj{njXL) ® (xljf^l) • (8) 

/=i 

Its non-locality can be expressed through a momentum-dependent formfactor J^{k), which is also a matrix in the 
Dirac s pace, attached to each fermion field. We will analyze the formfactors, including their dependence on density. 



in sect. IV B . After color-averaging, one obtains the effective interaction lagrangian 
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'4(7V2 

+ ^{^^^<y^.^r-T^f] , (9) 

where Nc is the number of colors and = (r, i) is an isospin matrix. In the pseudoscalar channel the interaction 
combines attraction for the isospin-1 (pion) channel with repulsion (due to the extra i) for isospin-0 (77'). Similarly, one 
finds attraction in the scalar isospin-0 (cr) channel (responsible for spontaneous chiral symmetry breaking) together 
with repulsion in the scalar isospin-1 channel (ao). 

In practice we will calculate correlation functions and the mean field effective potential in the Hartree-Fock approxi- 
mation. For this purpose, it is convenient to construct an effective s-channel kernel including the exchange term. This 
is made possible by the simple (separable) form of the momentum dependence. Using this kernel, one can reproduce 
the result of a Hartee-Fock calculation by evaluating the Hartree term only. In short-hand notation we will refer to 
the kernel as the effective meson or diquark lagrangian. From the Fierz identities given in appendix we obtain the 
following kernel for color singlet and octet qq states 



r — 



N, 



4(^2_i)(^-^'^"^M-^"-^^)'}' (10) 



■Z[I 

N, 



again being attractive in the a and tt channel, repulsive in the rj' and ag channel. Analogously, we can construct the 
effective interaction for color-antisymmetric 3 and -symmetric 6 diquarks. The result is 

= ^ {-^^ {(^^^T^CT^X\Ti,){^pT^X\CT*V 
^{'fT'^CT^X\^^T^){^^T^T2XYi,-,CT*^F)\ 
+ ^ {^^ T'^ CT^X%a^,T^,){^^T^r^X%a^,CT*^'^)^ , (11) 

where is the antisymmetric Pauli matrix, and Xa^s are the antisymmetric (color 3) and symmetric (color 6) color 
generators (normalized in an unconventional way, tr(A''A'') = Nc5'^^ , in order to facilitate the comparison between 
mesons and diquarks). In the color 3 channel, the interaction is attractive for scalar {ip'^Cj^il') diquarks, and repulsive 
for pseudoscalar (ip'^Cip) diquarks. 

As we have already discussed in our previous paper in the case of two-color (7Vc=2) QCD there exists an 
additional Pauli-Giirsey symmetry (PGSY) p3,p5| which mixes quarks with antiquarks. It also manifests itself in the 
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lagrangians given above, as in this case the couphng constants in qq and qq channels are identical, i.e., diquarks (the 
baryons of the iVc=2-theory) are degenerate with the corresponding mesons. Chiral symmetry breaking then implies 
that scalar diquarks are also Goldstone bosons, with their mass vanishing in the chiral limit {i.e., for current quark 
masses m = 0). 



B. 7-^-Molecule Induced Interactions 



Using the 't Hooft interaction introduced in the last section one can calculate correlation functions in a systematic 
expansion in multi-instanton interactions, starting from direct instantons graphs and proceeding to two-instanton or 
instanton-antiinstanton graphs, as well as more complicated clusters. In simple cases, like the set of RPA diagrams 



discussed in sect. II], one can sum a whole series of terms involving infinitely many instantons. But if the instanton 
ensemble is strongly correlated, this method may become very inefficient. In that case it is more useful to determine 
the effective vertex for a given cluster, and fix the strength of the vertex by calculating the concentration of clusters 
from the partition function. The simplest kind of cluster that can arise in the instanton ensemble are instanton- 
antiinstanton molecules. We have observed the formation of these clusters at high temperature and at large Nf 



in both analytic |^,^ and numerical simulations of the instanton ensemble |20 . In these cases, molecules are 
intimately connected with chiral symmetry restoration. An ensemble of molecules does not have delocalized zero 
modes or collective eigenstates, and the chiral condensate is zero. 

In the high density problem, the role of molecules is twofold. First, the concentration of instanton-antiinstanton 
molecules in the ensemble may be dynamically enhanced for similar reasons as in the case of high temperature or 



large number of flavors. We will discuss this problem in detail in sect. VIII. Second, the BCS instability is due to 
quark-quark scattering, or four-fermion operators. The 't Hooft vertex is a (2A^/)-fermion operator and does not 
automatically lead to an instability for Nf > 2. However, an instanton-antiinstanton molecule can always generate 
an effective four-fermion interaction, with the additional {2Nf — 4) fermion lines being internal. 

The effective four-fermion vertex induced by instanton-antiinstanton molecules was evaluated in . The result is 
particularly simple if the relative color orientation is fixed such that the instanton-antiinstanton interaction is most 
attractive. In that case one has 

CiA = G„,oi 1^ [(V^7pV')' + (^7m75^)'] - ^^^^^^ _ [(^7mA"^)' + (^7m75A'^V')'] 

[(^^7^^)' - (^7^75^-)'] - 2n!{N^- 1) [(^^^^»' - (^7^75 A»2] | . (12) 



Similar to the procedure leading to Eq. (^_1|), this interaction can be rearranged into an effective qq vertex. In the 
color antitriplet channel the result is 



oi [jfjj;^^) [{^'^Cj5T2X'X^){^^l5r2XW^) - (V'^Cr2A^V)(V5T2A^C^^)] 
+ ^^^(^^ _ ^) [{i''^Cj^l5r2yXi')i^Wr2\'XC^'^) - (V^C7^T2tA:^V)(V!7m^2tA1CV;^)] | . (13) 



Even in the case Nf — 2, there are two important differences as compared to the single-instanton vertex. First, 
since molecules are topologically neutral, the interaction is U{1)a invariant. This implies that it does not distinguish 
between scalar and pseudoscalar diquarks. Second, whereas the 't Hooft vertex only operates in scalar (and tensor) 
channels, molecules also provide an interaction in vector-meson and diquark channels. The coupling constant is related 
to the density of molecules and has to be determined from the partition function of the instanton liquid. We will 



study this problem in sect. VIII E 



III. DIQUARKS IN THE RANDOM PHASE APPROXIMATION 



Color superconductivity implies that the high density phase is composed of diquark Cooper pairs. In a weakly 
coupled BCS system, the expression 'Cooper pair' should not be taken too literally: it is not tightly bound and the 
range of its wave function is large compared to average inter-particle separations (i.e., the cube root of the inverse 
particle density, d = n~^/^). In QCD this is not necessarily the case. The gap can be quite large, and the existence 
of an intermediate phase of diquarks with or without Bose condensation is not a priori excluded. 
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For this reason we first study the possibility of diquark bound states in vacuum |2^^,^. In QCD with Nc > 2 
there are, of course, no gauge invariant diquark states. Instead, one can study correlation functions of heavy-light Qqq 
states, where the heavy quark Q serves to neutralize color. Effectively, this corresponds to a diquark correlator in the 
presence of a Wilson line. The "gauge invariant" diquark masses extracted from these correlators then measure the 
mass of the heavy Qqq state minus the mass of the heavy quark. In a dense medium, the Wilson line is not necessary, 
and color is neutralized by light quarks of the third color. 

In the effective fermionic theory derived in the previous section, diquarks can appear as physical bound states. 
These states should be interpreted as building blocks in the formation of baryons and dense matter. In practice, we 
study diquark correlation functions and look for poles in the diquark propagators at momenta \p\ < 2M, where M is 
the constituent quark mass. For simplicity, instead of the exact instanton formfactors ^{p), we employ in this section 
an euclidean 0(4) symmetric cutoff A, being adjusted to yield a realistic constituent quark mass M. 




FIG. 1. The quantity KJ entering the denominators of the T matrix for the scalar (S), pseudoscalar (P) and tensor (T) 
diquark channels. 



The Bethe-Salpeter equation for the two-body T matrix in a given channel can be written as 



T{q) ^ MC-'O^) X + ztr j 



F(p + q/2)(0'C)S^(-p + q/2)r(q) 



Introducing the notation 



k,k' 

Eq. (|lj) can be expressed schematically as 

T = K{1 + JT) . 

This equation has the solution 

T = (1 - KJ)-^K. 
In both diquark and mesonic channels J takes the form 



Jfc,fe,(p)=*tr j ^SF(p + q/2)0%(p-q/2)0'^ 



(14) 

(15) 

(16) 
(17) 

(18) 



since CSp{p)C ^ = Sf{—p)- From Eq. (^Tj) one has the three different operator structures, O = T2X'^i'-f5 (scalar 
channel), O = T2X\ (pseudoscalar channel) and O = r2A^cr^^ (tensor channel), which lead to 



Jss^-2h{M) + 2q^l2{q^M) 

JPP = 2/i(M) + 2(4M2 - q^)l2{q^,M) 

Jr^j. = -16M2/2(g^A^) 



(19) 
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with the two standard integrals 



They are readily evaluated for Euclidean momenta and analytically continued to Minkowski space to yield 

(22) 



27r2 



+ In ■ 



A2 + M2 



. 9 ^ ^^C 



M2 4M2-g2 g2 

m — IT --TT + 2\ / arctan 



A2 + M2 y (72 y 4^/2 _ ^2 

2A2 ^ / 4(A2 + Af2)-/ / g2 ^ 

4(A2 + Af2)-g2 j y ^5 ^^ct^'^ y 4(A2 + M2)-g2j ■ ^^-^^ 

The chiral gap equation in the scalar qq channel then becomes 

^ /i(M) = 1 , (24) 



8A^2 



which provides a relation between the coupling g and the constituent quark mass M. 

The conditions for the existence of poles in the corresponding channels of the diquark T matrix are 

(-2/1 (Af) + 2q^h{q^,M)) = 1 (25) 



87V2(^r^_l) 
for the scalar channel, 

- (^^ „i) (2-^1 W + 2(4M' - 9')/2(g^ M)) = 1, (26) 

for the pseudoscalar one and 

-^^ -(-16A/2/2(g^ A/)) = 1 (27) 

for the tensor one. If the l.h.s. crosses 1, there is a bound state. This is illustrated in Fig. |^ where we have taken 
M = 350 MeV, A = 900 MeV for definiteness. One finds that only the scalar diquark is bound, with a binding energy 
of about 200 MeV. 

This result agrees with numerical simulations of the instanton liquid |^^, which included all diagrams in this 
interaction. However, it is at variance with the conclusion drawn in ||25| , where no bound scalar diquark was found 
in the same model. We believe that the discrepancy is due to the fact that the authors of |Q only used part of the 
interaction in the diquark channel. In this work, we have performed a Hartree-Fock calculation with the full one- 
instanton interaction. A lattice calculation of diquark masses was performed in pl| . These authors find a significant 
scalar-vector diquark mass splitting, but no scalar diquark bound state. On the other hand, they also obtained a too 
small 7V-A mass splitting, and a mass ratio mj^/mp that is too large. 

One can also argue that there should be some continuity when going from the theory with A^c = 2 to A^c = 3. In the 
former case the scalar diquark is the partner of the pion and the vector diquark is the partner of p. This implies that 
the scalar diquark binding is large, Mdq,v — fndq,s = 'nip — m^r — 600Af eV . It is then natural that some remnant of 
the binding is left at — 3, since the coupling constant in the scalar diquark channel for Nc = 3 is only reduced by a 
factor of 2 as compared to the Nc = 2 case. This corresponds to a 'Wc = 2 + e" picture of the baryon octet: a tightly 
bound scalar diquark loosely coupled to the third quark. A number of phenomenological observations (reviewed, e.g., 
in j3^]) actually supports the validity of this picture for real QCD. The decuplet baryons, on the other hand, do not 
contain scalar diquarks, and therefore should be generic 3-body objects. This picture is quite contrary to another 
(and much better known) view of baryonic structure, the large Nc limit. Here both A'^, A as well as other members of 
the octet and decuplet are basically the same heavy object, slowly rotating with slightly different angular momenta. 
Indeed, as one reads off from the lagrangians given in the previous section, in this limit the diquark coupling tends 
to zero, and the scalar diquark binding disappears. 
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IV. INSTANTON-INDUCED INTERACTION IN DENSE MATTER 



A. Quark Zero Modes 

In the previous section the instanton-induced interactions between quarks were approximated by effective local 
4-ferniion vertices. In the microscopic treatment of sect. |l[ the external quarks couple to the quark zero modes in 
the instanton field, leading to a nonlocal profile function for the interaction vertex with a size characterized by the 
typical instanton radius of about /3=0.33 fm. In dense matter at zero temperature, the single-instanton solution itself 
is not affected by the surrounding quarks. The zero-mode wave functions, however, are density dependent leading 
to important modifications of the instanton-induced interactions in the medium. They can be constructed from the 
Dirac equation at finite (quark-) chemical potential, 

{i Pi - tfni)cf>i = . (28) 

The correct solution was obtained in p3| , p^ : 

A . P e^* / 2^ 2 + 7sin(H]e-^* 

^/ = «7; 5^- — 5 XL, 29) 

where the spinor xl arises from an antisymmetric (singlet) coupling of spin and color wave functions, as before. Note 
that the solution of the adjoint Dirac equation, 

</.|(a;;-//) (z^/-iAi74) =0 , (30) 

carries the chemical potential argument with opposite sign. This is necessary for a consistent definition of expectation 
values at finite /i, and in particular renders a finite norm. 



d X (i)\{x\-[i) (\)i(x;[i) = \ , (31) 

whereas without the extra sign one has 

d^x (/'|(a;; /i) 0/(2:; /i) = oo . (32) 



This singularity, corresponding to the particle-particle channel, is in fact directly related to well-known BCS singularity 
one encounters when resumming an effective (attractive) particle-particle interaction around the Fermi surface (see, 
e.g., ref. [||). 

B. Instanton Form Factors 

Using the explicit form of the zero mode wave function, we calculate the form factor from the Fourier transform 
^ = y" d4a;0(x)e~^'=-^ 

= 2ip \ dRB? \ drisin'^r) dOsine {-fodt + J ■ k cos Odr) 

Jo Jo Jo ^ 

(cos( ur) + - sm(ur))e~^^ -it,, a\ 

5~~ 5 ^ XL, [•J'i) 

-I- a;^ 

where x^ = r'^ and k = k/k. Introducing hyper-spherical coordinates for the integration, r = Rsinrj, t = Rcosrj, 
the result can be expressed through two scalar functions A{uj, k, fi) and B[uj, k, n), given in appendix H, as 

[loB{uj, k,p)+j- kA{u;, k, fi)]xL = ^Xl , (34) 
The finite density zero mode wave functions ^ have the symmetry properties 

ipiuj, k, p) — ip*{—u), —k, = ip*{uj, k, —p) . (35) 
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The combination J^{uj,k,iJ,) = ip*{Lij,k, ij,)Gq^{ui, k, fi) appears in the effective quark interaction on each propagator 
entering or exiting the instanton-induced vertex. In the mean field approach, when two of the propagators participating 
in the vertex have the same momentum, it is useful to combine them into two new formfactors. For a propagator 
entering the instanton vertex and another exiting with the same momentum one obtains 



a = J-{—uj, — fc, ii)''OJ-{uj, k, ^) , 
where O is a matrix with Dirac, color and flavor indices. For an overall unit matrix, one has 

a = - + il ■ k){'yaB*{uj, fc, ^) + 7 • kA*{uj, k, y)f 

(7oi(w - i^i) + ij ■ k) 
= {A*^ + B*^){k^ + {uj ~ ifif) 
= ar + ictj . 



(36) 



(37) 



For a propagator entering the instanton vertex and a transposed one exiting with the same momentum one finds 

13 = OT{-uj, -k, fifOT{uj, fc, fi) . (38) 
When the Dirac part of O is C75, where C is the charge conjugating matrix, 

P = Clzhoii-i^ - im) - • k){'JoB*{-uj, k, ^) - 7^ • kA*{~uj, k, ^)) 

Cj5{loB*{uj, k,fi) +j ■ kA*{uj, k, ^))(7oi(u; - i^J.) + ■ k) 
= (7o(w + i/^) + 7- k){'yoB{u!,k,fi) + 7- kA{uj, k, ^j.)) 

{'yoB*{uj,k,fi) + 7 • kA*{uj,k,fi))(jQ{uj -ifi)+j-k) 
= (cj2 + k^ + /i^)(|^P + \B\^) + 2fiki{A*B - AB*) 



+i-foi ■ k[2fik{\A\^ + \B\^) + (u^ ■ 
= Pr + ilol ■ kPi . 

We note that a and [3 have the same symmetry as in Eq. 



k' + ^i')i{A* B - AB*)] 



(39) 





FIG. 2. The moduli squared of the form factors, \a\^ = a'f. + (top panels) and |/3p = f3f. — j3f (bottom panels), for two 
difterent values of the chemical potential: /i = (left column) and ^ = 300 MeV (right column) . 
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C. Quark Overlap Matrix Elements 



The instanton fomifactors discussed above are designed for momentum space calculations, in particular for extract- 
ing effective interactions between quarks in the mean-field f rame work. However, in the statistical mechanics treatment 
of the instanton liquid partition function presented in sect. VIII, the coordinate space description is the more suitable 
one. For that we will need the explicit form of the fermionic overlap matrix element Tja which at finite density takes 
the form 



Tia{z,u; h) 



4{x 



zi; -y) [iip - ifi'jA:) (t)A{x - za; 
- zi; ~n) (i ^~ (t>A{x - za; m) 



(40) 



The second line is again obtained by virtue of the Dirac equation when choosing the sum ansatz for the gauge- 
field configurations, Afj, = Aj^ + A^. Tja plays a crucial role in the fermionic determinant of the instanton partition 
function, where it generates the fermionic interaction ('quark hopping amplitude') between Ps and A's and is therefore 
responsible for correlations in the instanton liquid. In particular, Tia controls the probability of forming molecules. 

The definite chirality of the zero modes (in the limit of vanishing current quark masses) entails that /-/ and 
A-A matrix elements are zero. In the vacuum Lorentz invariance implies that the overlap matrix element can be 
characterized by a single scalar function |]36|| , e.g., Tja = i u ■ z f{z). In the medium this is no longer true and Tja 
must be calculated in terms of two independent scalar functions /i, /2 according to 



Tja{z,u\ijl) 



i U4 /i(T,r;^) +i /2(T,r;/i) 



(41) 



They are shown in Fig. ^, see also 38|. Similar to the finite temperature case, we observe a strong enhancement 
with increasing n in the temporal direction. Moreover, the fermionic interaction becomes very long range, ^ [i^ jx^ (at 
finite temperature it was limited to the Matsubara box of size 1 jT enforced by periodic boundary conditions) . In the 
spatial direction, the exponential damping exp[— TrTr] in the finite-T case is replaced by oscillations ~ sin(/ir). The 
latter also effectively suppress the hopping amplitude once the r-integration in the partition function is performed. 
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FIG. 3. Quark-induced /-^-interaction at finite density for the most attractive color orientation Ui = l, u=0 as well as r—Q 
(left panel) and for U4=0, = l and r=0 (right panel). 



From the strong enhancement of Tja with increasing fi one may already anticipate the relevance of mo lecule 
configurations at finite densities ||3^]. This issue will be quantitatively investigated in the 'cocktail' model in sect. VIII . 



V. THE TWO-FLAVOR PROBLEM 



This section will be devoted to study finite density two-flavor QCD in th e mean-field approximation using the exact 
momentum dependent instanton profile functions as discussed in sect. IV B As density increases the basic competition 
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will be between the chiral condensate {qq) and the scalar ud diquark condensate in the {qq) channel, representing the 
color supercond uctin g state as discussed in Refs. We will first discuss the corresponding coupled gap equations 



formalism (sect. V A) and then proceed to the numerical results and their interpretation in sect. VB 



A. Mean-Field Grand Canonical Potential at Finite /i 

For the evaluation of the grand canonical potential we employ the Cornwall- Jackiw-Tomboulis (CJT) effective 
action, which is elucidated in more detail in Appendix^. It involves 3 type of propagators (including their antiparticle 
pendants) corresponding to single (anti-) quarks carrying color charge that participates in the diquark condensate 
(Gi, Gi), single (anti-) quarks carrying color charge that is not part of the diquark condensate (G2, G2), and (anti-) 
diquarks {F F). The minimization of the action with respect to (w.r.t.) these propagators generates the following six 
gap equations: 



-{Gi-FGi^F)-^ + G^^- 


- Mia 





-{Gi-FGi^Fy^ +Gq^ + 


Mia* 


= 


— G2 ^ + Gq ^ — 


M2a* 


= 


— G2 + Gq 4- 


M2a* 


= 


Gi^F{Gi-FGi^F)~^ 


+ iAp 


= 


Gi^F{Gi~ FGi^F)-^ H 


-iAp* 


= 0, 



(42) 



where 



Ml = 2g 
M2 = 2g 
A^2g 



1 



1 



(tr(Gi + G2)a) 



1 



-(trA8(Gi + G2)a) , 



8iV; 



V2 16N2(N2 - 1) 
2 (tr(Gi + G2)a) - -4 .«^!!m2^ ,^ (trA8(Gi + G2)a) 



1 



8N^{N, - 1) 



V2 16N2(N2-1)' 
tr(FC75A2r2/3), 



(43) 



are the two chiral masses and the diquark gap, and Gq (Gq) is the bare (anti-) quark propagator defined through 
Eq. (C2). As before, the traces involve momentum integrations. The chiral masses are linear combinations of the qq 
condensates, tr(Gia) and tr(G2a), while the gap A is proportional to the qq condensate, tr(FC75A2r2/3). Eqs. ( |4^ ) 
represent a coupled system of gap equations in the chiral and diquark masses, Mi,M2 and A. To determine their 
solutions, one needs to know the explicit form of the propagators Gi,G2 and F. They are constructed from the 
coupled set of Eqs. (^2|). Using the relation G{p) = —G^{~p) and the transposition property G757J — 'fp.Gj^, one 
can rearrange Eqs. (p2[) into Gorkov-type equations (note that G,F,F do not commute) as 



G2{p) = Go(p) 
Gi{p)^Go{p) 
F{p) = F{p)Mia 



G{p)M2a{p)Go{p) 
G{p)Miaip)Go{p) 



F(p)(zAG75A2T2/3(p)Go(p) 



{p)G^{-p) + G{p)iAI3*{p)G^,\2T2Gl{-p) . 



(44) 



A graphic representation of these equations is displayed in Fig. IJ. 
The Gorkov equations can be solved in algebraic form yielding 



G2{p) ^ {G,\p) ~ M2a{p))-^ 

Gi{p) - (Go i(p) - Mia{p) + A''p*{p){G^\-p) - Mia* {p))-^ (3{p)y 
F{p) = iAGi{p)p*{p){G^\-p) - Mia*{p)r^G-„\2T2 
= iA{G^\p) - Mia{p)r^l3*{p)Gi{-p)G-i5\2T2 

= F(p)G"f^\2T2 

F(p) = iAG^.,X2T2{G^\-p) - Mia*ip))-'Pip)Giip) 
= iAGT.\2T2Gii-p)P{p){Go'{p) - Mia{p))-^ 

= G^z^2T2F{p) ■ 



(45) 
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G^P) 



Go(p) 



Go(p) 
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F(p) 
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Go(p) F(p) 



5 "■a "2 
+ 



Go(-P) Gi(p) Go(-p) F(p) 

FIG. 4. Diagrammatic representation of the Gorkov Eqs. (^^. 

To obtain the thermodynamic state variables such as pressure and energy density we need to know the explicit 
dependence of the thermodynamic potential on the mass p aram eters. This can be achieved by reinserting the explicit 
solutions of the propagators into the effective action, Eq. (C12), constituing the grand canonical potential times the 
4-volume, —V4^l{Mi, AI2, A). For that purpose we evaluate the propagators more explicitly. After some algebra, one 
can rewrite Gi from Eq. (Uq) as 



Gi(p) = l^° 



a 



flavor 



~«7o 



[lo - i^i)((uj + tfif + fc2 + Afi*2) 



-Ml 



((w + inf + fc2 + Mf) + M*A''iPl - (3f)] } 



(46) 



where Mi{p) — Mia{p), and 

V 



|(^ _ j^)2 + fc2 ^ ^^2)|2 ^ ^4(^2 _ ^2)2 _ 8k^^A^'PrP^ 
+ 2A2(cj2 + fc2 + /i2)(/32 + + 2|Afi|2A2(/32 „ ^2)^ 



(47) 



Using the relation (C5) for the Dirac part (tro) of the trace-log in the kinetic part of Eq. ( C12| ) (here the trace does 
not include the momentum integration), one finds 



-tro ln(-GiGi + GiFGj"^F) = -41nX», 



(48) 



and from Eq. (| 



Re [trD(Go ^Gi - 1)] 



92\2 



Re[((a; + ifiY + k')Mt] + |Afi|* + A^(/3,^ - pf) 



4/cAiA^/3,A + A\lo' + k' + fi^)if]t + K) 



Re [trD(Gia)] = - — 



+2|Mi|2a2(/32_/32 
4 

Ml 



Re[((cj + ^^Jif + e)Ml] + |Afi|4 



Re 



trD(F/3) 



4 
A 



A\0', - 0if - Akp^^'Mr + A'(^' + k'+^i') 



(49) 

(50) 
(51) 
(52) 



The analogous quantities involving G2 and M2 are obtained from the above by substituting A — ^ 0, Afi —f M2. 
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Using the above relations our final expression for becomes 



fi(Mi,M2,A) 



duj- 



Pdk 



27r3 

4Re [trD(Go^Gi - 1)] + 8 



4 In I? - 2 In I - ifif + + M^)\^ 



|(tJ-i/i)2+fc2+M|)|2 

A{y^ d^__Re[trD(2Gia + G2a)]} -^{y^ rf^__Re [trD(Gia - G^a)] } 

trD(F/3)] }' . (53) 



du „ '^^ Re 



27r3 



The global minimum of f2 w.r.t. Mi, M2, A at each /i determines the thermodynamically stable phase and the values 
of Ml , M2 , A are the chiral masses and color superconducting gap in that phase. The extrema of ft at each /i can be 
found by equating the derivatives w.r.t. Mi, M2, A to zero. Equivalently, one can verify that after differentiating the 
expression for w.r.t. Mi,M2 and A one recovers Eqs. (43). Their solutions correspond to the local extrema of 
and represent different possible phases. We shall discuss them in the next subsection. Phase transitions correspond 
to two distinct minima of Q having an equal value (first order), or merging together (second order). 

One should recall that in this fomulation the coupling constant g is an integration variable: an inverse Laplace 
transformation was used to exponentiate the instanton vertex. However, in the thermodynamic limit the saddle-point 
approximation for then- integration becomes exact (since it is multiplied by the 4- volume V4). Identifying the potential 
energy of fl in Eq. (M) as —gU, the integral in question is 



Z (X dg exp[—gU 



y Hg)] , 

where N/V is the total instanton density. The saddle point is then found to be at 

_ N 1 

9'max Yj ' 



(54) 



(55) 



VIII within a 



Thus, at the saddle point the new potential energy is —N/V\n{U), up to a constant. 

The real question is how to determine the /i dependence of N/V. This will be addressed in sect, 
statistical mechanics treatment taking into account correlations in the instanton ensemble. It will be shown there that 
the simplifying assumption of a constant total instanton density is indeed reasonably justified. Another approximation 
concerns the density-dependence of the second key parameter, the average instanton radius p, which defines the scale 
in all instanton calculations. Lacking better knowledge, we also assume that it does not vary significantly at the 
chemical potentials under consideration. 



B. Nf = 2 Phase Diagram 

Solutions of the ga p e quations are extrema of fi, but only minima represent a thermodynamic phase. There are 4 
types of solutions to (|4^): 

1. There is a chiral condensate, but no diquark one, i.e.. Mi — M2 — M, A = 0, 

2. There is a diquark condensate, but no chiral one, i.e.. Mi — M2 = 0, A ^ 0, 

3. Both condensates are present, i.e., all Mi, Af2, A 7^ 0, 

4. No condensates at all, i.e.. Mi = AI2 = A = 0. 

The last case corresponds to free quarks and it is easy to see that, at T = 0, having at least one condensate is always 
more favorable. The rather complicated expressions for the formfactors require the integrals in Eq. (^3|), together with 
Eq. (|5^), to be evaluated numerically for each value of fi varying from to 500 MeV. We shall call the three types 
of solutions described above phase 1, 2 and 3, and fix our two input numbers as p = 1/3 fm and {N/V)p^ = 0.0116, 
corresponding to N/V = 0.94 fm^**. 

With these values we find that at fi — the system is in phase 1 with the familiar value of about 330 MeV for 
the chiral mass M at /i = 0. We normalize the grand-canonical potential such that the minimum at /i = has zero 



13 



pressure, f2 = 0. We find minima for phase 1 in the range of fi from to 360 MeV, for phase 2 for all values of /i and 
for phase 3 only for values of fi between 250 and 290 MeV. The values of for all three phases are shown in Fig. ||. 





phase 1 




phase 2 - 




V ^ phase 3 _ 




\ 

\ 

s 







phase 1 
phase 2 
phase 3 




200 300 
fj. (MeV) 



270 
M (MeV) 



FIG. 5. Q, for phases 1, 2 and 3. The right panel is a magnification of the region between 250 and 290 MeV, where phase 3 
exists. 



Phase 1 dominates until 250 MeV where the system makes a transition to phase 3, and then, at 288 MeV, there 
is a transition to phase 2. The coupling constant g for all three phases changes little (not shown). The mass M for 
phase 1, the gap A for phase 2 and the two masses Mi, M2 and the gap A for phase 3 are shown in Fig. O. 




FIG. 6. The left panel shows M for phase 1 and A for phase 2. The right panel is a magnification of the region between 250 
and 290 MeV, with the masses and gaps of all three phases. 



To understand the nature of the phase transitions we study the dependence of CI on the three parameters Mi , M2 
and A. As seen in the right panel of Fig. ^, phase 3 starts at the values of the masses and the gap of phase 1 
at fj, — 250 MeV. This is an indication of a second order phase transition, as the corresponding chiral and diquark 
condensates (which are proportional to Mi, M2 and A) are the first derivatives of O w.r.t. these masses and the gap. 
The emergence of the phase transition at /i = 250Mey s exhibited even more clearly in the left panel of Fig. 
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FIG. 7. The left panel shows profiles of SI at = 270 and 246 MeV along the direction that at ^ — 270 MeV connects the 
solution for phase 1 = 0) and phase 3 = 1). The right panel shows f2 at /i = 284, 288 and 292 MeV along the direction 
that a.t jj, — 288 MeV connects the solution for phase 1 = 0) and phase 3 — 1). 

At n = 270 MeV, we have plotted the values of on a straight line in the parameter space (Mi = Mi^phi + 
^6Mi,M2 = M2,phi + ^SM2,A = Aphi + ^SA), connecting phase 1 = 0) with phase 3 = 1). We see that the 
solution for phase 1 is not a minimum but a local maximum of fi. The second curve shows ft just before the transition, 
where phase 3 has not yet emerged, but the minimum is quite flat indicating that the second derivative of Q is close 
to 0. When moving along the same direction {SMi,SM2,6A) in the parameter space we see a classic example of 
a second order phase transition, when at certain value of the parameter /i the symmetry changes (the SU(3) color 
symmetry is broken to SU(2)), a diquark condensate appears, the second derivative of Q goes through and the old 
solution (with the higher symmetry) turns from a minimum into a maximum. The new absolute minimum is phase 
3. Of course any new extrema can only turn up in pairs, but all solutions are symmetric w.r.t. the sign of A so that 
the two new minima must be at ±A. 

The second phase transition at /i = 286 MeV is analyzed in the right panel of Fig. |^. The middle curve shows the 
values of at /Lt = 270MeV on a straight line in the parameter space (Mi = Mi^ph2+£,5Mi^ M2 = M2^ph2+^5M2, A = 
Aph2 + ^SA), connecting phase 2 (^ = 0) with phase 3 (^ = 1). At this value of /i phase 2 already dominates over 
phase 3. The opposite situation is observed from the upper curve which is for /i = 284 MeV with the cross-section of 

in the same direction in the parameter space. Obviously a first order phase transition occurs for some fj, between 
these two values. However, the barrier between the phases is quite low and the values of the parameters (Mi, M2, A) 
(and the condensates) of the two phases are quite close, so it is a rather weak first order transition. This is further 
supported by the fact that at /i = 290 MeV the minimum for phase 3 disappears quite rapidly (by going through an 
inflection point), as seen from the right panel of Fig. ^. For the lowest curve only phase 2 exists, but the remnant of 
the inflection point is visible. 

There are some shortcomings in the mean-field analysis as presented here. Below some critical /Xc, which marks the 
onset transition, no physical quantity should depend on fj,. Due to the fact that the instanton zero modes explicitly 
depend on fi - however small - this is not respected in our calculation. Nevertheless, the variation of M and Q below 
the phase transition is quite small. The result might be further improved by taking into account the dependence of the 
instanton density and size on /i. Another problem is that the onset transition happens quite early, at /j, ~ 250 MeV, 
whereas the expected onset corresponds to a third of the nucleon mass minus the binding energy of nuclear matter 
(939 — 16)/3MeV ~ 308 MeV. Again, this might be related to various approximations employed. 

It is interesting to note that we do find all three phases to exist, not just the chirally broken and superconducting 
phases, but also a phase with chiral symmetry breaking and diquark condensation. The existence of the latter phase is 
not a very firm prediction as the difference in energy density of the three phases in the transition regions is rather small. 
In fact, phase 3 was not observed in the NJL calculation of |^ or the instanton calculation of |Q or ||l3|. The latter 
work uses slightly different techniques to evaluate the grand canonical potential. We will return to the (spe culative) 



phase with simultaneous diquark condensation and chiral symmetry breaking in the discussion in sect. IX B 
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VI. THREE FLAVOR QCD IN THE CHIRAL LIMIT 



The situation becomes more involved if one includes the strange quark. Since the critical chemical potential 
He ^ 300 — 350 MeV is larger than the strange quark mass iris — 140 MeV, strange quarks have to be included 
whenever there is time for strangeness to equilibrate. There are several qualitatively new features in going from 
Nf = 2 to Nf = 3. First, since Nf = Nc, there are new order parameters in which the color and flavor orientation of 
the condensate is locked Second, the instanton-induced interaction is a six-fermion vertex, so it does not directly 
induce the BCS instability. 

In this section we will consider Nf — 3 flavor QCD in the chiral limit. In that case, we expect that in the low density 
phase chiral symmetry is broken by a quark condensate (uu) — (dd) = (ss). In the high density phase, quark pairs are 
condensed. One possible form of ordering is the analog of the Nf = 2 diquark condensate, {qfC^^q^) — A^eijk£abc- 
This order parameter breaks color SU{3)c ^ SU{2)c and the chiral 5[/(3)l x SU{3)r 5L/(2)l x SU{2)r. A more 
attractive possibility is provided by the following order parameter |^ 

iirRlj%) = liCh^PRr" {AlS^aSjb + A2S^bS,a) ■ (56) 

Here, Pr is the projector on right-handed quark fields, and there is an analogous expression for left-handed fields 
also. This order parameters breaks color and chiral symmetry down to the diagonal subgroup SU(3)c+l+r- Since 
the color symmetry is completely broken, there is a gap in the spectrum for all 9 quarks and 8 gluons. This already 
suggests that the phase characterized by Eq. (^6|) should be preferred over the Nf = 2 like phase, in which only 4 
quark states are gaped. We will see this more explicitly in the next section. 

The order parameter ( |5^ ) breaks chiral symmetry since the residual symmetry couples flavor rotations of right and 
left handed quarks. The diagonal symmetry acts on the quark fields as 

qf ^ {U*)^,U^y„ (57) 

where U is an element of SU(3). The most general form of the quark condensate that is consistent with this symmetry 
is 



= 2 iP^r y [^^ ~ 3 ^« ) ^^'^ + 2S8<5",<5';. j . (58) 

At zero density we expect Eg to be zero, but in the high density phase both Sq ^^nd Sg will in general be non-zero. 
It is important to note that even though the above argument establishes that chiral symmetry is broken, it does not 



show how a chiral condensate is actually formed. From the superfluid order parameter (56), we can directly form the 
chiral order parameter ((^l^a)^) ~ {qhqh) {qRqR) , but not the chiral condensate (qLqR)- This is because ( |5^ ) violates 
right (and left) handed quark number by two units, whereas (qLqR) violates right and left handed quark number 
by one unit. In other words, the order parameter leaves a discrete chiral symmetry unbroken, and this symmetry 
prevents the quark condensate from acquiring an expectation value. But this discrete symmetry is explicitly broken by 
instantons. In the color-flavor-locked phase, we can saturate four of the external legs of the instanton vertex (qLqR)^ 
with the condensate and obtain an effective interaction (qLqL) (qRqR) (qLqR.) which leads to the formation of a quark 
condensate. 

In the case of three massless flavors, the 't Hooft interaction is a flavor antisymmetric six-fermion interaction |4 

1 / 2iVc + 1 

GN.iN^-lf^'^'^'''''''' \2N, + 'i^ 

3 



£ = G&{2'Kpf ^ -'^hf2h<^gig2g; ( , ^ (^^./i V'fl.gi )(V'L,/2 V'B,g2 )(V'L,/3 V'fl.ga) (59) 



lPB.,g,) + {L ^ R) . 



In the following, we will consider the somewhat more general case of a U{1)a violating six-fermion interaction char- 
acterized by two independent coupling constants Gg i and Gg 2 , corresponding to the scalar and tensor terms in 
Eq. (H). _ 

In the vicinity of the Fermi surface, six-fermion terms are suppressed w.r.tfour-fermion interactions. This is the 
Cooper phenomenon: Near the Fermi surface, the only interaction that is not suppressed is 2 ^ 2 scattering, where 
the two particles are back-to-back. In the more modern language of the renormalization group one finds that the 
strength of the six-fermion interaction is reduced as one integrates out states away from the Fermi surface [|4|-Q . 
In the context of the mean-field approximation employed here, we will see that the gap equation has a logarithmic 
enhancement in the case of four-fermion interactions, but not for six- (or even higher) fermion vertices. 



16 



For this reason we will have to consider the effect of four-fermion interactions. We already stressed that instanton- 
antiinstanton pairs provide a four-fermion interaction for any number of flavors. In terms of left- and right-handed 
fermions, the interaction is 

-^{^LlMii'B.lM - ^[^^2 ~_ 1^) (V^lTmA^V'l) (Vii?.7M A>;,J I . (60) 

In the following, we shall study the condensates ( p6| ) and ( p8| ) for an interaction given by the sum of the four-fermion 
vertex ( |60| ) and the six-fermion vertex (|59|). In this section, we will consider the coupling constants G4 and G g to be 



arbitrary parameters, constrained mainly by the known value of the quark condensate at zero density. In sect. VIII E 
we shall try to determine these couplings from the partition function of the instanton liquid. 

The system of gap equations for the three-flavor case can be derived along the same lines as the two-flavor case 
discussed in the previous section. However, the resulting equations are algebraically much more involved. In order 
to keep the presentation reasonably simple, we will ignore the instanton form-factors and take the interaction to be 
point-like. As we saw in the last section in the case of Nf = 2, this approximation does not qualitatively affect the 
results. 



0.1 



> 

< 



0.01 



0.001 




FIG. 8. Chiral and superconducting gaps ai,2 and (5i,2 as a function of the chemical potential for the three-flavor model in 
the chiral limit. 



We shall calculate the thermodynamic potential in the mean-field approximation. This calculation is somewhat 
complicated by the fact that the color-flavor structure of the propagators is quite involved. The first step is to 
determine the quadratic part of the action in the mean-field approximation. For this purpose, we close off all except 
two legs of the interaction. The result is 



G4 ( y So - ^J^s) + G6,i Ua^I + 8E0S8 - y E^ ^ 9^2^ 



G6,2 144E^ + 48EoE8 - 168E^ - 30Aj 



G4^E8 4- Ge,! ( -24EoE8 + 10^1 - -A^ 



G, 



6,2 



144EoE8 -I- 120E| -t- 18A^ 



'-Gi + Ge.i (3So - 2E8) + Ge,2l2 (3Eo - 4E8) 



(61) 
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We note that the interaction is only sensitive to the antisymmetric part A^i = Ai — A2 of the {qq) order parameter. 
This is different from the OGE interaction considered in |^], but does not make much of a difference in practice, 
since even in that case the solution of the gap equation has As/A^i <C 1, where A5 = Ai + A2. We also note 
that there is a chiral symmetry breaking qq interaction proportional to A^. This is as expected: Color-flavor-locking 
combined with instantons leads to chiral symmetry breaking. We also stress that both ingredients, instantons and 
color-flavor-locking, are essential. 

The {qq) and {qq) mass terms can be diagonalized simultaneously. In general, the mass term can be decomposed as 

+ .91M1} + (g5^,,C754,j) {hMi + /2M2} , (62) 

where we introduced the color-flavor matrices 

Afo-<5"''*y, Mi=5'Jj''^-, Ah^S'^^S^. (63) 

The matrices AIi, M2, commute. This means that there is a color-flavor basis in which ( |62|) becomes diagonal. 
We denote the quark fields in this basis by cjjp, with p = 1, . . . , 9. The mass term becomes 

y , {(Tl {4>px4>pm) + h {4>pM.C^b4>p,R)] I + {(72 {(1>9,l4>9,r) + (^2 (09,flC7509,i?)} , (64) 
\p=l ) 

where ai = go, 6i = —f2 is eightfold degenerate and cr2 = go + Sgi, 62 = 3/i + /2- In our case 
<Ji = G4 (^y Eo - l^s) + Ge.i ("84^2 + SEqEs - y E^ + ^A^ 



Gea + 48EoE8 - 168E| - SOA^^ (65) 



y Eo + y Eg^ + Ge.i (^84E^ - 64EoE8 + y5 

+ Ge,2 (144E2 - 384EoE8 + 192E2 -|- 24A^) (66) 
61 = ^62 - A^ 0G4 + G6,i (3Eo - Eg) -I- G6,2 (-3Eo + 4E8)^ . (67) 

The potential for the mean field is given by closing of all external legs of the interaction. This way we get two-loop 
graphs proportional to G4 and three-loop graphs proportional to Gg. In the mean- field approximation, we have 



F = G4 2A3i + 12E2 + -I]2 + Ge.i 504Eg - 384EoE| + — E| + 24 (3Eo - 2E8) A\ 



320. 

+ G6,2 (864Eg - 2304EoE| + 1280E| - 144 (3Eo - 4E8) A^ j . (68) 

In the quadratic part of the interaction we can now integrate over the fermion fields. Since the color-flavor structure 
is already diagonal, we get a sum of nine terms, each (in principle) with different gap parameters. We finally obtain 
the following result for the free energy 

F = -8e{ai,Si)-e{cr2,S2) + V . (69) 
Here, the single particle energy is given by 



'^^'^^^ / (2^ \ \/iEp-f^r+S' + sJiEp + nr+S^ \ , (70) 



(2^)3 

and Ep = + a"^. The mean- field parameters Eo,E8, A^ are determined by making the free energy stationary 
(9F)/((9E,;) = {dF)/{dAi) = 0. This gives three coupled gap equations that have to be solved numerically. 

Before we present the results we have to discuss how to fix the coupling constants G4 and Gg. We take G1.2 to 
have the relative size implied by the instanton interaction (^9|). If we were to ignore random instantons, and only 
had instanton-antiinstanton pairs, the four-fermion interaction would break chiral symmetry for G4 > 7.5 A^^. We 
consider this to be the upper limit on this interaction. In order to see how large the gaps in the three-flavor case can 



18 



possibly be, we take G4 just below this limit G4 = 7.4A~^. Gg is then fixed by the requirement that for nis = 150 
MeV (see next section) we get a reasonable constituent u, d mass of 400 MeV. This gives Gq = 12. OA^^. 

Results for the various gaps are shown in Fig. |[ Note that the superconducting gap is smaller as compared to the 
two-flavor case. This is because diquark condensation is now due to pairs, not individual instantons, and we restricted 
the size of the corresponding coupling such that it does not lead to chiral condensation at /i = 0. This is similar to 
the scenario of Alford et al. where superconductivity for Nf = 3 is driven by one-gluon exchange. Again, reason 
dictates that the corresponding coupling is below the critical coupling for chiral condensation. 

Instantons lead to chiral condensation in the diquark condensed phase. This is immediately clear because if we 
take the six-fermion vertex and close off four legs by two diquark insertions, the remaining interaction violates chiral 
symmetry. Color-flavor locking is nevertheless essential. For the two-flavor superfluid order parameter, instantons 
only lead to a non-zero (ss). Alford et al. realized that chiral symmetry would be broken, but could not calculate 
the size of the effect in their model. Here we flnd it to be very small. The maximum constituent mass generated 
in the diquark condensed phase is less than 10 MeV. Qualitatively, this is not hard to understand: the constituent 
mass arises from terms in ( |6l| ) that are proportional to A^. These terms arise as exchange terms from the original 
interaction (jsol), so they are suppressed by degeneracy factors 2NfNc. In addition to that, the constituent mass 
is driven by the superconducting gap squared, which is already about an order of magnitude smaller than the zero 
density chiral gap. 

There is one more important direct instanton effect in the high density phase. For Nf — 3 all chirally invariant 
four-fermion interactions (molecules, OGE, etc.) are U{1)a invariant, and do not distinguish between scalar and 
pseudoscalar diquarks. This means that a parity broken vacuum characterized by the order parameter 

{qtCq'j) = ^^6^a5bJ + A2S,bS,a (71) 

is degenerate with the parity conserving vacuum considered here. The same is true for an arbitrary linear combination 
of positive and negative parity condensates. The degeneracy is lifted by the six-fermion interaction in conjunction with 
flnite quark masses or non-vanishing chiral condensates. This implies that the difference in energy density between 
the parity broken and parity conserving vacua is small. This is different from the Nf = 2 case, where the four-fermion 
interaction distinguishes between { fn\ ) and ( p6|), and the energy difference is big. This effect is also different from the 
scenario considered by Pisarski and Rischke |]47[, who argued that the parity broken vacuum is degenerate with the 
parity conserving one if instanton effects are small. In three-flavor QCD in the chiral limit parity broken and parity 
conserving vacua are almost degenerate, even if instanton effects are not small. 



VII. FLAVOR SYMMETRY BREAKING 



The situation is even more complicated if we take flavor symmetry breaking into account. For simplicity, we will 
restrict ourselves to m„ = rriii = and rrig 7^ 0. It is clear that as — > 00, we have to recover the two-flavor scenario, 
with the order parameter given by 

{q^Cj^q'j) - Ande^.se'^'' • (72) 

Note that in the two-flavor case the color orientation of the condensate is arbitrary, but for three flavors the choice ( [7^ ) 
is preferred because it preserves an SU{2) subgroup of the diagonal SU{3)c+l+r.- We might also consider additional 
gap parameters that have a different color orientation, but the corresponding gap equation simply decouples and the 
solution (except in the limit ttt-s — > 00) is not energetically favored. 

Since flavor symmetry is broken, the structure of the quark condensate is also more complicated. The following 
ansatz generalizes Eq. (pq): 

i&tj) = ^iPRT" ((so - ^Ss) 5""% + S.<5'"'<5.3<5,3 + 2J:sS''^S) + Sg^iFi + Eg^sPs) , (73) 

where Pi = S'^^SaS^'^Sjs and P2 = d'^^SaS) + S'\S''^Sj3- 

There are a number of complications that occur once flavor symmetry is broken, and it is hard to take into account 
all of these effects at the same time. In the following we will concentrate on the dynamical interplay between a flavor 
symmetric four-fermion interaction generated by one-gluon exchange or instanton pairs and the flavor symmetry 
breaking four-fermion vertex that comes from the six-fermion 't Hooft interaction and a strange mass insertion. In 
addition to that, we have to take into account that there is no pairing between strange and non-strange quarks if the 
mismatch between the Fermi momenta is too big. The BCS instability arises for pairs with total momentum zero 
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where both of the individual momenta are on the Fermi surface. This is not possible if the masses are different and 
the Fermi surfaces are shifted. In the presence of pairing the Fermi surface is not sharp, but smeared out over an 
energy range given by the gap. This means that pairing between strange and non-strange quarks is suppressed if the 
mismatch between the Fermi momenta exceeds the gap, 



ml/{ApF) > A 



(74) 



Moreover, there is the problem that the color-flavor matrix characterizing the most general diquark mass term does not 
commute with the color-flavor structure of the (qq) mass term. This means that we cannot simultaneously diagonalize 
the two mass terms, and write the free energy in the simple form (^3|). Instead it seems unavoidable to deal with the 
full (spin, color, flavor, and (qq) versus (qq)) matrix structure of the quark propagator. On the other hand, we found 
that, except possibly in a small regime, quark and diquark condensates do not coexist below the critical chemical 
potential, and that the quark condensate in the high density phase is small. In the following we will therefore treat 
the quark condensate in the high density phase as a small perturbation. 

The important new ingredient if 7^ is the presence of a four-fermion interaction which operates exclusively in 
the u,d quark sector. This interaction arises from the 't Hooft interaction, Eq. (|59|), by closing off two external legs 
by a strange mass insertion. The result is 



1 



2N,{N^ - 1) 



(75) 



which (of course) has the form of the Nf — 2 't Hooft interaction, but with a coupling constant controlled by the 
parameter 3ms/{4:TT'^p^). So, unlike in the OGE-based works, the value of the strange quark mass has not just 
kinematical but also dynamical significance. 

Our model then consists of a flavor symmetric four-fermion interaction, the flavor symmetry breaking four-fermion 
interaction (|75|), and the flavor symmetric six-fcrmion interaction (|59|). In order to compare with work of ARW [Q, 
we take the flavor symmetric four-fermion interaction to be one-gluon exchange. We could equally well have used the 
instanton-antiinstanton induced interaction - qualitatively this makes very little difference. 

In this model, the mass term becomes 
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XSo + (7G4,i + 6G4,2) So + (84G6,i -I- 144G6,2) {^l + SqS. 
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K^, - (7G4,i -I- 6G4,2) So - (84G6,i + 144G6,2) SqS, 



K 



{2Aa - As) 



^'0 



-j(2Aa + As) 



^3ab 



-KAud 



1 



(G4,i - 12G4,2) {Aa + 2Aud) 



3 2 

The quark mass term is already diagonal, while the diquark mass term is of the form 

K.fil^oq'hJ {/iMi + /2A/2 + /3M3} 
with Mo. 1,2 as before and M3 = e^'^^e^ij. This mass matrix has four eigenvalues, 

■51 = ±/2 
82 = ±(/2 - h) 
1 



33,4 



3/1 + 2/2 + /3 ± V9/? + 2/1/3 + /; 



(76) 



(77) 



(78) 
(79) 

(80) 



with degeneracies di = 4, 3, 1, 1. The free energy is given by F = — di£{ai, Si) + V as before, where the potential is 
V = 16K (3^2 + 2SoS, + + 6 (7G4,i + 604,2) Sq + (504G6,i + 864G6,2) (So + S, 



4if 3Ai-3A|-f4A2, 



AAaAuo) + (g4,i - 12G4,2) [Aa + 2A, 



(81) 
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An important difference as compared to the flavor symmetric case is that we have to take into account the kinematic 
restriction discussed above. In channels that involve pairing between strange and non-strange quarks we restrict the 
integration to the regime {Ep — fi)^ > (rn^/{App) — 6'^). For a more detailed discussion, we refer the reader to |4|] and 
p9| . Again, the gap equation follows from the requirement that F is stationary w.r.t. the gap parameters and Aj. 




FIG. 9. Superfluid order parameters and gaps as a function of ma for /i — 350 MeV (left panel) and /i — 500 MeV (right 
panel). The upper panels show the up-down and up-strange components A^d and Aus of the color-flavor locked state. The 
lower panels show the maximum, minimum, and average superconducting gap. 



Numerical results are shown in Fig. ^. The parameters were fixed as described in the last section. In the figure we 
plot Aud and Aus as a function of rus for two different chemical potentials /i = 0.35 and 0.5 GeV. In particular we 
show the value of the largest, the smallest, as well as the average gap. In the two-flavor case, the smallest gap is zero 
and the average gap is 4/9 of the maximum gap. In the three-flavor case, the smallest gap is 1/2 and the average gap 
5/9 of the maximum gap. 

We observe that there is a sharp transition between the two-flavor scenario (A^s = 0, Aud 0) and the three-flavor 
scenario {Aud = A„s ^ 0) that takes place around the physical value of the strange quark mass, m"'** = 65 MeV for 
/i = 0.35 GeV, and mf"^^ = 160 MeV for fi — 0.5 GeV. The ratio Aud/ ^us grows roughly linearly already for small 
TOg. This is different from the results of and an instanton effect. As discussed above, instantons induce a 

four-fermion interaction among light quarks which is proportional to m^.. We should note that we have not included 
the possibility of a dynamically generated contribution to the strange quark mass in the superfluid phase. In terms 
of the current mass, this effect will shift the critical mass to smaller values. 



VIII. STATISTICAL MECHANICS OF THE INSTANTON LIQUID 



A. The Cocktail Model at non-zero Density 



In this section we take a step beyond the mean-field approximation in that we allow for possible clustering effects 
in the instanton ensemble. This is achieved by employing a somewhat different formalism. In the previous sections we 
started from an effective quark interaction obtained by integrating out the gluonic (instanton-) fields in the underlying 
partition function. In this section we reverse the strategy and integrate over the fermion (quark) fields first. This 
leads to the following partition function for the instanton ensemble at finite density: 

Z^nst{t^) = J2 N \NJ. n / "(z^^) P^' n 'ietl^ P + imj- z/i74) . (82) 

N+,N- +' ■ 1=1 •' /=1 
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The original QCD path-integral over all possible gluonic field configurations has been converted into an integration 
over the collective coordinates flj = {zi, pi,ui\ (position, size and color orientation) of A'^_|_ instantons and 
antiinstantons. The single-instanton amplitude n(/9/) contains the semi-classical tunneling rate (including one-loop 
quantum corrections), as well as the Jacobian arising from the introduction of collective coordinates. The instanton 
interactions can be divided into a gluonic part Sint and a fermionic part represented by the determinant of the Dirac 
operator. It is usually approximated in the subspace of zero modes, i.e., 



det(« p — i^J,j4) 



dot ( T,A^^) ^ 



(83) 



where Tja is the fermionic overlap matrix element discussed in sect. [VC, see Eq. (^0|). When restricted to the zero 
mode basis, the fermionic determinant is in fact equivalent to the sum of all closed loop diagrams to all orders in the 
't Hooft effective interaction. The non-hermiticity of the finitc-/i Dirac operator is reflected by the fact that 



TA/(/i) = T;^(-/i) 



(84) 



i.e., the fermionic determinant is complex, entailing the well-known 'sign'-problem in the partition function, which 
will be addressed below. 

In the statistical mechanics treatment chosen here, spontaneous chiral symmetry breaking in the vacuum is generated 
by randomly distributed uncorrelated anti- /instantons which allow for a delocalization of the associated quark quasi- 
zero modes corresponding to the formation of a nonzero (qq) condensate state. In other words, quarks can travel 
arbitrarily long distances by randomly jumping from one instanton to another (antiinstanton) , and thus may carry 
their chiral charge to spatial infinity where it effectively becomes 'lost'. 

In this picture chiral restoration can in principle proceed in two ways: either instantons disappear altogether, or 
they rearrange into some finite clusters which do no longer support any finite (qq) condensate. In the limit of very 
large temperature or density instantons will eventually disappear, because Debye screening of the large gluon fields 
inside the instantons leads to a strong suppression of the tunneling rate. Nevertheless, in the case of finite temperature 
QCD it was argued |2^] that this cannot be the relevant mechanism for chiral restoration, since lattice simulations 
have observed the transition at rather low temperatures of ~ 150 MeV. This, in turn, led to the suggestion that 
chiral restoration proceeds through the formation of I- A molecules. Further support for this idea is provided by lattice 
measurements of the instanton density at finite temperatures, showing no depletion below and a smooth onset of 
the expected Debye-screening above . More recent studies ^1 50 1 seem to find more direct indications for I- A 

molecule formation at T Tc, but their quantitative role in the transition remains to be clarified. 

At finite density additional possibilities for clustering of the (chirally asymmetric) random instanton liquid into 
chirally symmetric configurations are available. In particular, random instantons can be supported by diquark con- 
densates for Nf = 2 or by a combination of diquark and quark-antiquark ones for larger Nf. However, as we will show 
in this section, effects of I- A molecule formation may still be an important element in the finite-/i chiral restoration 
transition, as was shown in Ref. In some sense this is not really surprising: when both T and /i are sufficiently 
large so that the all condensates are absent {i.e., in the true QGP phase), I-A molecules should constitute the preferred 
configuration for any remaining instanton component in the system. 

To investigate the interplay between the various components in the finite density partition function more quantita- 
tively we resort to the 'cocktail-model' introduced in Refs. [^,^. Here, the instanton ensemble is decomposed into 
a mixture of random ("atomic") and "molecular" configurations, which yields a grand canonical partition function of 
the form 



inst 



Na,Nrr 



Nal NJ. 



N„ 



(85) 



In the thermodynamic limit V4 ^ oo, and using the Stirling formula, the free energy (thermodynamic potential) 
becomes 



Qa-t-m 
inst 



ln[2^°+'"l 





- Tim. In 















(86) 



The atomic and molecular 'activities' are i5ll,Ef: 
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= C ^^("-4) 1 [T,^(A.)rA/(/.)]^^ P^^^ . (87) 

The underlying approximation in this approach is that the values of the hopping amplitudes Tja in each individual 
configuration are replaced by a product of their mean square values in an uncorrelated ensemble. To establish a 
connection to the (chiral) quark condensate and constituent mass, one can use the mean-field approximation to 
express them via the "atomic" density na as 

1 /3 



im) = -- [-^n. ) (88) 



M = -Cm I {TTpfiqq) = -Cm {ttp) ^^n, , (89) 

respectively. In the original work of Ref. |^^, where these relations have been first derived, the coefficient Cm = 1, 
leading to an effective quark mass of M* ~ 200 MeV. The latter is to be understood as the average mass of a constituent 
quark at finite (euclidean) momentum participating in tunneling processes associated with instantons. With increasing 
four-momentum M* is appreciably reduced and therefore does not correspond to the usual constituent quark mass 
M (defined at zero momentum). We account for this by using Cm — 2, yielding M — 400 MeV. The minimization 
of the total f2 over Ua and rim determines their equilibrium values in the ensemble for given T, p. The corresponding 
gap equation for the constituent quark mass is the direct analog of the mean-field equation conjugated to the (qq) 
condensate, but expressed in terms of different variables. 

For a refined treatment at finite densities we supplement the cocktail model by two additional components. Following 
the arguments given above, we have to account for the possibility that the random instanton component can become 
engaged in diquark chains, first observed in numerical simulations of the instanton ensemble in the high-density 
limit [p8[ . The pertinent term in the free energy reads 



Vld{nd\p) = -n^ln 



ezd 



nd 



(90) 



with the associated activity 



Zd = 2C p"-^ e-^-* f ^ j d^z du [TjA{p)TjA{p)h P^j ■ (91) 

The subscript "3" indicates the color projection in some predefined direction characterizing the color vector of the 
diquark. In the same way that Ua determines the constituent quark mass M, the superconducting gap is related to 
the density Ud as 



In analogy to Eq. (^) we use Ca = 2, but also perform calculations with Ca — 1-5 to assess the inherent uncertainty 
of this mean-field estimate. Notice that in contrary to Eqs. (^^, the fermionic overlap matrix element enters Eq. ( pl| ) 
as (Tja)'^, which, in fact, causes the z-integration to diverge. This is precisely the BCS-singularity of an attractive 
interaction in the particle-particle channel, here encountered in coordinate space. The standard procedure to treat this 
singularity is to start from a new ground state which a priori has the gap built into the fermion propagators, thereby 
regulating the integrals. The net effect of the gap on the overlap integrals is a damping factor for the intermediate 
quark propagators, which is delineated in appendix 

The second refinement consists of including a Fermi sphere of quark-'quasi-particles' in the free energy, representing 
the contribution of quark non-zero modes. The final expression for the thermodynamic potential then becomes 

nina, nm, Ud; p) = ^l+T+'^ina, n^^Ud; p) + f^^^^^^jAf, A; p) , (93) 
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where the quark contribution is simply given by 



^ZarkiM, A; ^i) = e,(Af, A; /i) - n,(M, A; ^l) 



(94) 



with 

eq{M; = gq J d^k ujk npifJ- - Wfc) , (95) 

gq = 2NcNf and = Af^ + (in the superconducting phase the dispersion relations of paired quarks picks up an 
additional dependence on the gap A, see below). 

At this point it is instructive to compare the 'cocktail model' with other approaches. The (simpler) philosophy em- 
ployed in the first part of this paper introduces the multi-fermion interactions with a constant coupling g, independent 
of condensates and temperature/density. However, since g is itself generated through instantons, density variations of 
the latter will affect the coupling. A step towards including this feature was made by Carter and Diakonov [|l3|: their 
coupling parameter, called A (which essentially corresponds to Za without the fermionic determinant part), was sub- 
jected to minimization. The resulting gap equation relates the mean instanton density Ua to A and the condensates. 
However, instead of calculating A and then finding (as done here), they stayed in the mean-field approximation in 
which the instanton density remains unchanged. Our results to be discussed below support the (approximate) validity 
of this assumption. 

To investigate possible mechanisms for chiral symmetry restoration at finite density in some detail, in particular the 
competition between diquark and molecule formation, we will in the following separately discuss various versions of 
the cocktail with increasing complexity, i.e., the two-flavor model including (anti-) instantons, molecules and a Fermi 
sphere of constituent quarks (sect. VHI B| ), additionally incl uding t he simplest ud-pairing as discussed in Refs. 
(sect. VHI C ) , and the three-flavor case (sect. VHI D ) . In sect. VIII E we also give estimates for the density-dependence 
of effective coupling constants for molecule-induced (anti-)quark-quark interactions, which naturally emerge from the 
formalism employed in this section. 



B. Two Flavor Cocktail Model without Diquark Condensates 

In this subsection we basically follow the approach of Rcf . , generalizing it to flnite density . The issue here 
is to assess the potential role of I- A molecule formation in chiral symmetry restoration at finite density, without the 
additional complication of superconducting gaps. 

For the actual calculations we now have to face the problem of the complex fermionic determinant appearing in the 
various activities, Eqs. ( p7| ) and (|9l|). As has been suggested in Ref. |^^, it can be solved under the assumption that 
the gluonic interaction does not exhibit a pronounced dependence on the color angles, approximating it by an average 
value (see below). As a result, the color dependence in the activities only enters through the combinations of T/yi's, 
which then can be integrated analytically, rendering the fermionic determinant real. For two flavors one obtains 

Za{z4,r) CX / duTiA{fJ.)Tj^{~fl) = ^[/i+Jf + /2^/2"] 



z„(z4,r)cx / rfu[T/A(/i)r;^(-A*)] ' = TTTn^ (96) 



^'^ " AN,iN^-l) 

where = fi{zL^) are defined through Eqs. ( ^ ) and (^). Note that the color averaging of the former complex 
expressions is sufficient to yield real- valued activities. Also note that, whereas ZmizA, r) is a positive definite quantity, 
this is not the case for Za{z4, r) due to the oscillations in fi^2{r). In fact, when further integrating Za{z4, r) over space- 
time (the remaining four collective coordinates), delicate cancellations occur, which require accurate numerical values 
for the otherwise one easily encounters negative/incorrect results for both activities at finite chemical potential. 
The gluonic interaction entering Eqs. ( |87| ) has been approximated by an average repulsion Sint = Kp'^{na + 2n,„), 
where k = (3 /2p^{N /V), j3 = b/2 + 3Nf/A — 2, 6 = — |-/V/- The free parameter k characterizes the diluteness of 
the ensemble. In the case Nc — 3 and Nf — 2 we have chosen k ~ 130 in order to reproduce the phenomenological 
value for the diluteness of the instanton vacuum. The normalization constant C oc (Aq^/))^ can also be fixed in 
the vacuum by requiring that the absolute minimum of the thermodynamic potential flinstip- = 0;na,rim) appears 
at a total instanton density of N/V ~ Ua + 2n„i = 1.4 fm^^, being realized for na=1.34 fm^'* and n„j=0.03 fm~^. 
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which gives Aqcd — 260 MeV. The smallness of the molecular component in the vacuum is a consequence of the 
large entropy associated with quantum fluctuations of the color angles, randomizing the system. Available lattice 
data agree that such correlations are indeed small at T = /i = 0. 
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FIG. 10. Results for the two-fiavor cocktail model without diquark pairing; the left part shows the densities of random 
instantons ('atomic' component) and I-A molecules (upper left panel), the pressure p — —Q, (middle left panel) and the 
constituent quark mass (lower left panel) after minimizing the free energy, Eq. (DSl), w.r.t. Ua and Um. In the right panel the 
free energy is displayed as a function of the constituent quark mass M oc nj , indicating a first order transition from the 
minimum at finite M to the one at M = 0. 



Our numerical findings for the Nc = 3, Nf — 2 atomic+molecular instanton cocktail model at finite density are 
summarized in Fig. ^ At small /i essentially nothing happens until, at a critical value /ic — 310 MeV, the system 
jumps into the chirally restored phase, the latter being characterized by = 0. The transition is of first order, as 
can be seen by inspection of the Af -dependence of the free energy (right panel of Fig. [lO| ). Below the pressure 
actually decreases slightly with increasing fi indicating a mixed phase-type instability, similar to what has been 
discussed in Refs. 53|. The total instanton density at the transition (residing in /-A- molecules) is appreciable, 
N/V = 2n,„ ~ 1.1 frn"^^, providing the major part of the pressure at this point. In other words: a substantial part 
of the non-perturbative vacuum pressure persists in the chirally restored phase. The vacuum pressure of p(/i = 0) = 
0.6 GeV fm~'^ (indicated by the dotted line in the middle left panel of Fig. M) is not recovered until a chemical 
potential of /zq — 350 MeV, corresponding to a free quark number density of ri^ = 1.15 fm~^ (naively, this translates 
into a nucleon density of n^r = 2.4no with the normal nuclear matter density tlq — 0.16 fm~^). For all chemical 
potentials in between, < fj, < fj,Q, the system is mechanically instable, possibly indicating droplet formation in the 
region M < ^ < fiQ (as suggested in Ref. Q), where, with a finite quark density, the pressure p(/x) is below its vacuum 
value. 

A further comment concerning the numerical value of the critical chemical potential, which is very close to a third 
of the nucleon mass, is in order. Given the various approximations applied it should be regarded as a coincidence. At 
the same time it most likely provides a lower bound for the true value, as has been the case for the finite temperature 



25 



calculations of Ref. |Q (resulting in a critical temperature which is about 20% lower than observed on the lattice and 
might well be related to the fact that the instanton model lacks explicit confinement). On the other hand, the inclusion 
of superconducting gaps may further reduce the critical chemical potential and in fact dominate the transition, as will 
be discussed in the following section. 



C. Two Flavor Cocktail Model Including Pairing 

We now address the effects of diquark condensation. For clarity, let us first ignore the I- A molecule component 
in the cocktail to study the competition between chiral and diquark condensates only. The additional inclusion of 
a pairing gap resulting from ud diquark formation at the Fermi surface as discussed in Ref. Q modifies the quark 
quasiparticle contribution according to 

r!,^(M)=trlog[D(M,A)] . (97) 

Here, D denotes the quasiparticle quark-propagator now including the BCS gap A. All interaction contributions are 
effectively accounted for through the instanton part ^f^gt^ ^1^. (|90|), (^l|). The resulting expression for the free 
energy then becomes 

niua, na- = 17^+?(n„, n^; + ^ [2f^f (M, A; + {N, - 2)17^^(M; ^)] (98) 

(the last term accounting for unpaired quarks), which now has to be minimized w.r.t. and nj^. The results for two 
different values of the vacuum instanton density and the (not precisely determined) coefficient Ca for calculating the 
pairing gap (cf. Eq. (p2[)) are summarized in Tab. | We find that color superconductivity appears at critical chemical 
potentials around /Zc — 300 MeV, very similar to the values of the previous section where only I- A molecule formation 
was considered. The associated gaps range between 120-180 MeV. These results are consistent within 20% with the 
findings of sect. ^ and those of Ref. We should also note that the calculated gaps in our original work ||] are 

significantly smaller (below ^ 100 MeV) as those were effectively obtained for 2-1-1 flavors, i.e., in Ref. we included 
the effect of a reduced (constituent) strange quark mass Mg in the closed-off strange quark loop of the six-fermion 
instanton vertex, which decreases the effective instanton-induced coupling constant for the four-quark interaction 
by about 60% in the chirally restored phase. Another feature that emerges here is that the total instanton density 
changes little across the transition, which a posteriori justifies to assume it as constant in the mean-field calculations 
of sect. 

As in the previous section, there is an intermediate constituent quark-diquark phase, similar to what was found in 
sect. but here it again has small negative pressure which makes it mechanically unstable against the formation of 
a mixed phase. 

Let us now turn to the full two-flavor cocktail model with simultaneous account for the chiral and diquark conden- 
sates as well as I-A molecules. The total free energy 



r!(n,,nd,n™;^) = "(n,,, n^, n™; /^) + ^ [217^ (M, A;^) + (A^, - 2)r!^^(A/; /i)] , (99) 

limized w.r.t. Ha, rim and rid. The results are shown in Fij 
(which most closely resembles the results of the calculations in sect. 



is to be minimized w.r.t. ria, rim and rid. The results are shown in Fig. [T l j, using the coefficient of Ca — 1.5 in Eq. (|9 





Ca 


fic [MeV] 


MO [MeV] 


A(Aie) [MeV] 


1 


2 


300 


355 


158 


1 


1.5 


265 


290 


120 


1.4 


2 


310 


375 


188 


1.4 


1.5 


270 


305 


142 



TABLE I. Parameter dependence of the critical chemical potential for chiral restoration in a cocktail model with chiral and 



diquark condensates (no I-A molecules included). 
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FIG. 11. The full two-flavor cocktail model including I-A molecules and diquark pairing; upper panel: densities of random in- 
stantons ('atomic' component, full line), I-A molecules (dashed line) and instantons engaged in 'diquark chains' (dashed-dotted 
line); middle panel: pressure p(/i) = — 0(/i), maximized over na,nm and Ud at each value of fj,; lower panel: constituent quark 
mass (full line, using Cm = 2 in Eq. (|9|) to give M(/i = 0) = 400 MeV) and diquark gap (dashed-dotted line, using Ca ~ 1-5 
in Eq. (^), as discussed in the text). 

As to be expected from the previous analysis, there is a delicate competition for chiral restoration between random 
instantons engaged in diquarks and I-A molecule formation. The former do, in fact, induce the chiral transition for 
all parameter ranges considered. 

Towards higher densities it may happen that a second transition occurs within the chirally restored phase char- 
acterized by a substantial jump in the molecule density and an accompanied drop of the diquark gap, (its location 
is somewhat sensitive to parameter choices). The singularity in Zd, on the other hand, guarantees that there is al- 
ways a finite (qq) condensate present, albeit possibly strongly reduced in molecule-dominated phases. In Fig. |l^, the 
transition occurs at /ic — 270 MeV into a chirally broken diquark phase, which, again, is mechanically unstable. The 
vacuum pressure is only recovered at fiQ = 295 MeV, i.e., the combined effect of diquarks and molecules further lowers 
the ^Q-values found when including only molecules (Fig. ^ where /xq = 350 MeV) or only diquarks (Tab. |, where 
Ho = 305 MeV). 
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Another important point to note is that the total instanton density, N/V = + + 2nm, is indeed essentially 



independent of the chemical potential within ~15% or so. E.g., in Fig. 10, the upward jump in the density of molecules 
is close to half the drop in density of the random instanton component. A similar continuity is seen in Fig. |l| for 
the more realistic case of the transition between the two kinds of random instanton liquids, associated with the two 
condensates {qq) and {qq). After all this is not too surprising since, at least for chemical potentials /i < 0.4 GeV, the 
by far dominant fraction of the total energy in the system is carried by the gluonic component residing in instantons, 
and the free energy itself must of course be continuous across the transition. 

Finally, a comment concerning the finite temperature behavior is in order. On the zero-density and finite-T axis 
it has been shown that molecules drive the chiral restoration. This implies that at finite density, starting from low 
T, the role of molecules should become more and more relevant; at the same time, color-superconducting gaps are 
well-known to be suppressed. Thus we expect that for any value of /i where there is diquark condensed phase for 
T = 0, a rearrangement into a molecule-dominated phase should occur when raising the temperature. As indicated 
by our zero temperature results, for chemical potentials /i > 300 MeV or so, the corresponding Tc(/i) line in the phase 
diagram might actually reach down to fairly low temperatures. 

D. Chiral Restoration for more Flavors 

As was already mentioned in the introduction, in vacuum the tendency towards chiral restoration in the instanton 
model strongly increases with the number of fiav ors, le ading to a chirally symmetric vacuum state for Nf as low as ~ 5. 



The reason for that has been discussed in sect. II B: the increased number of quark lines enhances the interactions 
between instantons and antiinstantons, making the random liquid less favorable. More specifically, the integral over 
color orientations within the molecule, being proportional to ((cos^)^''^-'), increases more strongly than that for the 
averaged 'random'-instanton configurations raised to the iV/-th power, {(cos6Y)^f. The former integral is strongly 
peaked at = 0, creating a "locking" of the color orientation within a molecule. 

From continuity one may thus expect that the critical chemical potential for chiral restoration should be further 
reduced when moving from two to more massless flavors. Here we would like to pursue the question in how far a third 
massless flavor impacts the results of the two-flavor case. For simplicity, we now ignore color superconductivity and 
consider an interplay between random and molecular components only. As discovered in |^ and further elaborated 
in the first part our article, starting from Nf = 3 a color-flavor locking phenomenon sets in, leading to a complicated 
set of {qq) and {qq) condensates. However, they are relatively small in magnitude, and can therefore be neglected for 
the present purpose. 

Compared to the Nf=2 case, the color integration for Nf—3 is substantially more involved. Using the appropriate 
relations for the integration over a string of six SU(3)-color matrices (see, e.g., Ref. Isl]), one obtains 



i{z4,r) cx / du[TiA{fJ.)Tj^{-fi)]^f 



16iVc(7Vc + 2)(iV2 



- (3iv, m)Hfr? + {f^nf2r] + [4 - n,] [1^/2 fi^') ■ im 



Note that Zm{z4,r) in the three-flavor case has no definite sign before integration over space-time (similar to the 
Nf=l case, i.e., Za{zi,r) from Eq. (p6|)), which inevitably entails partial cancellations. Apparently, only for an even 
number of flavors the space-time integrand of the finite density activities is positive definite. Performing again the 
minimization procedure in Ua and rim for the thermodynamic potential we find that the critical chemical potential is 
indeed further reduced, by about 10% to /Uc — 270 MeV (as compared to 310 MeV for Nf=2). This indeed complies 
with the above mentioned expectation that at a sufficiently large number of flavors a purely "molecular" vacuum 
is more preferable than a "random" one, and thus chiral symmetry would be unbroken even in vacuum. As this 
phenomenon reflects itself also along the /i-axis, we expect fic to be further reduced at Nf=A, possibly crossing zero 
at Nf^b. 

E. Molecule-Induced Effective Couplings 

In this section we apply the cocktail model to a microscopic estimate of the density-dependence in the effec tive 
coupling constants for (anti-)quark-quark interactions. In the mean-field framework employed in sects. VII 
such density-dependencies were not accounted for. 
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To begin with let us recall the expression for the effective coupling constant for single- instanton induced interactions 
as, e.g., used in Ref. 



G^nst = / dpTi{p,n) p^f {2npY 



with 



i{p,p)=Cn^ {Sn^/g^f exp[-87r^/g{pf] p- 



(101) 



(102) 



denoting the single-instanton distribution, i.e., the semi-classical tunneling amplitude. Replacing the p-integration 
by an average value for j?, one has 



G.nst = n± {pY'f (27rp) 



(103) 



with the anti-/instanton density n^. Notice that (for Nf = 2) six additional powers of p appear in the integral of 
Eq. (101) as compared to the usual expression for the density, — J dp n{p). This entails a significantly larger 



average value p than the typical instanton size of about p ~ 1/3 fm. In order to be consistent the standard value 
of the zero-density constituent quark mass M ~ 400 MeV (which requires Ginst — 19.1 fm^), one should have the 
effective size saturating this integral to be p = 0.51 fm. 

The effective coupling constant for molecule- induced interactions has been first derived in Ref. p8| , where it is 
written as 



Gmoi = / n{pi,p2) dpi dp- 



2 (2^Pi)' {2np2 f 
with the total molecule amplitude 

n{pi,P2) = du d^z n{pi) n{p2) TiA{u,zf^i pf^ p^' 



(104) 



(105) 



which is nothing but the molecular activity given in Eq. (p^). The graphical interpretation of Eq. (104) is quite 
transparent: starting from the molecule amplitude, where all 2Nf quark legs are closed within the molecule, one 
'cuts' open two of them (corresponding to the division by Tfj^) which provides an effective interaction between four 
external (anti-) quarks. Inserting (105) into (104), and replacing again the size integrals by using average values for 
p we obtain 



Gri 



{2t: pfp^^h 



2Aff-2 



du d z Tia{u, z) f 



(106) 



We see that for Nf — 2 the individual p integ rations carry additional powers of p**, while it is p^ for Nf — 3. Those 



powers should be compared to p^ in Eq. (101) and p° for the usual instanton density. As these last two integrals lead 
to p ~ 0.51 and 1/3 fm, respectively, we interpolate between them (lacking a mor e accurate determination of the size 



distribution), i.e., estimate the appropriate average size values entering Eq. (106) to be p ~ 0.43 fm for Nf = 2 and 



p ~ 0.47 fm for Nf — 3. The corresponding zero-density values for the cou plin g constants are (including a factor of 
16 accounting for the difference in color coefficients between Eqs. ( |l0| ) and (p^): 



G^nstip = 0,p - 0.51fm) = 19.1 fm^ 
16 G^i'^^ip = 0,p = 0.47fm) = 0.86 fm^ 
16 G^^^^{p = 0, p = 0.43fm) = 0.044 fm^ 



(107) 



Obviously, G„.st ^ 20 X (16 G'^^J^), and (16 G^if) ~ 20 x (16 G;„„, 
averaged over positions and sizes of an /-A-pair, is given by ||l5|| 



27r2 Np-^ _ 1 
iN'r~V~ ~ 40 



'). Recalling that the fermionic matrix element. 



(108) 



(for p — 1/3 fm), one readily understands the decreasing magnitudes of the coupling constants in terms of the 
additional powers of dimensionless diluteness of the ensemble entering Eq. (106). Note that the effect of color 
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integrations favors molecules, and the corresponding factors compensate, to some degree, powers of the small diluteness 
parameter. 

We are now in position to evaluate the density-dependence of the effective coupling constants. U sing th e fcrmio nic 



VIII D we 



matrix elements at finite /i, Eg. ([40| ), and performing the color integrations as discussed in sects. VIII B , 
obtain the results shown in Fig. |12|; we recall that throughout this paper possible effects from the Debye-screening of 
instantons at high densities have been ignored. Therefore Ginst is in fact a constant as it does not depend on T/^(/i). 
On the other hand, the behavior of the molecule-induced couplings depends on the number of flavors: whereas in the 
Nf = 2 case it decreases with fi, the opposite is found for Nf = 3. This is directly related to the /x-dependence of the 



activities calculated in sect. VIII B , since (/i) cx Tp^ (corresponding to Za) and G^^ (/i) oc T/^ (corresponding 

to zZ'~'^)- 

It is instructive to compare the values of the coupling constants with the perturbative OGE interaction. For 
large-angle scattering the coupling is 

GOGE = , (109) 

where Peff is some effective momentum transfer averaged over the Fermi sphere. With Peff — 0.5 — 1 GeV and 
as{peff) ~ 0.3 we have Gqge = 0.15 — 0.7 fm^. This is comparable to the effect of molecules for Nf = 2, and 
significantly larger than that for Nj = 3 at low /x. One should note, however, that the structure of the interactions is 
different. For example, OGE is flavor independent whereas instanton-induced interactions are flavor antisymmetric. 
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FIG. 12. Instanton-/molecule-induced effective 4-quark coupling constants as a function of chemical potential. 

IX. MORE PHASES, OUTLOOK AND EXPERIMENTAL CONSEQUENCES 

A. The phase diagram 

In this section we would like to put the main results obtained in this work into perspective and discuss the emerging 
picture of the QCD phase diagram. 

The main feature of the phase diagram for two-flavor QCD in the chiral limit and at T = is that at a criti- 
cal chemical potential /ic the system undergoes a transition from the chirally broken phase to the superconducting 



^For small-angle scattering there appears an additional logarithmic enhancement, which becomes relevant for color supercon- 
ductivity at asymptotically high 
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phase ||i|,p,|T^. With instanton- induced formfactors we also found a small window in which a chirally broken 'con- 
stituent quark' plus diquark condensate phase may exist. Unfortunately, this is not a robust prediction of t he mo del, 
since the free energy differences between the phases are small (in fact, in the cocktail model analysis of sect. VIII this 
phase is mechanically unstable). We will discuss this issue from a slightly different perspective in the next section. 

In any case, the mean-field approach predicts a strong first order phase transition, either from the vacuum phase 
to superfluid quark matter phase, or from the chirally broken diquark phase to the superfluid phase. This implies the 
existence of an inhomogenous phase at intermediate density, with dense quark matter bubbles immersed in the chirally 
broken phase. Of course, a more refined treatment should reproduce the fact that matter clusters into nucleons and 
nuclei. In the high density phase chiral symmetry is restored, but color-5J7(3) is broken to SU{2). We have not 
explored the possibility of further breaking SU{2) via color-6 condensates. These condensates seem to generate very 
small gaps for the up and down quarks of the third color ^ . 

For three massless flavors we also flnd a first order phase transition from the chirally broken vacuum phase to the 
superconducting phase. In spite of the fact that instanton-induced dynamics are very different from OGE considered 
in Ref. we also find that the preferred order parameter in the superconducting phase exhibits color-flavor- locking. 
Both color-5[/(3) and chiral SU{3)l x SU{3)r are broken, while the diagonal SU{3)c+l+r is preserved. But even 
though color-flavor locking in general implies that chiral symmetry is broken, instantons are crucial for generating a 
non-zero value of (qq). In practice, the value of the chiral condensate turns out to be small. 

Furthermore, we have considered the more general case of QCD with 2-f 1 flavors allowing for a massive strange 
quark. Of course, the two cases Nf = 2,3 discussed above emerge as limiting cases for rrig — > and —^ oo. For 
nis — 2\/A/i pairing between light and strange quarks becomes impossible, and there is a phase transition between 
the color-flavor locked phase and the two-flavor superconductor. Again, a signiflcant difference between the instanton 
model and schematic interactions abstracted from one-gluon exchange appears. In the case of OGE, the strange quark 
mass has a purely kinematical effect. Instantons induce four-fermion interactions of the type ms{ud){ud), which can 
generate large asymmetries between the the {ud) and (its) = {ds) components of the color-flavor locked state even for 
TOs < rrf"^ . Similar effects are well known from hadronic spectroscopy. 
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FIG. 13. Schematic QCD phase diagram in the chemical potential /i - temperature T plane. The small T-/^-region corre- 
sponds to ordinary hadronic matter, with broken chiral symmetry. The point M (from "multifragmentation" ) is the endpoint 
of the nuclear liquid-gas phase transition. The point E indicates where the first order line either terminates in a second order 
endpoint (for mu,Tnd ~ 0) or disappears (for finite light quark masses). CSC2 and CSC3 label the Nf = 2 and Nf = 3-type 
superconducting phases. The hypothetical intermediate quark-diquark phase is indicated by QDQ. 

In the present work we have not addressed the effects of finite temperature. One would expect that the first order 
chiral phase transition at /i ^ will persist for some range in temperature, until the transition becomes second 
order at a tricritical point In BCS theory one can also estimate that superfluidity disappears at a critical 



temperature Tc(/i) ~ 0.6A(/i,T = 0) (however, as discussed in sect. VIII the impact of instanton-antiinstanton 
molecules presumably reduces the BCS-coefficient appreciably). As explained above, the boundary between the two 
types of superconductors should be approximately determined by the condition ([T^). Since, asymptotically, the gap 
is expected to slowly grow as a function of chemical potential, the critical temperature will also grow. As a function 
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of chemical potential the system will eventually reach the color-flavor locked state for any value of the strange quark 
mass. On the other hand, as a function of temperature (at given /x) one might expect that, since the gap in the 
two-flavor superfluid is somewhat larger than in the color-flavor locked phase, the system first makes a transition to 
the two-flavor superconductor, and then to the quark-gluon plasma phase. Combinig these conjectures leads to the 
schematic phase diagram dsiplayed in Fig. |l3| (where we have also indicated the nuclear liquid-gas transition line). 



B. More Phases? 



Finally, we would like to make a few comments on the range of applicability and the limitations of our approach. 
First, we have restricted ourselves to an instanton model. Even though many of our conclusions apply to any 
phenomenologically successful effective interaction, many of our numbers are indeed model-dependent. 

Furthermore, we have made extensive use of the mean- field approximation (MFA). The MFA is valid provided 
the condensates are sufhciently smooth, and their fluctuations can be neglected. This assumption becomes better 
at large density, because in this limit the coupling is weak and Cooper pairs are large compared to the inter-quark 
distances (as is the case in ordinary superconductors). In the opposite limit of small density the simplest phase one 
can generate in the mean-field approach is a chirally asymmetric Fermi gas of constituent quarks with masses on the 
order of 400 MeV. Remarkably enough, our approach predicts that this phase is unstable against separation into a 
mixed phase, with high density quark droplets separated by pure vacuum (5^,01 . While this result appears to be very 
suggestive, the MFA cannot predict the actual composition of the clusters. 

Of course, we know that the "correct" clusters are nuclcons. Below nuclear matter saturation density, nucleons 
themselves will form a mixed phase of clusters (nuclei), but for larger density one has homogeneous nuclear matter. 
In the vicinity of nuclear saturation density [un ^ tiq ~ 0.16 fm^^), the equation of state is known experimentally. 
The behavior of the energy density as a function of density is commonly parameterized as 

e = (mjv - (5mAr) njv + ^^7^ ( 1 - — ) , (110) 

1» \ Uq J 

where SniN ~ 16 MeV is the binding energy (per nucleon) of nuclear matter and the compression modulus is on 
the order of K^^ — 200 — 300 MeV. This means that nuclear matter is rather stiff, that is, the pressure grows very 
fast as a function of density. The physical reason for the steep rise in pressure is not just Fermi motion, but, more 
importantly, a strong repulsive core in the nucleon-nucleon interaction. The standard lore is that the steep rise in 
the pressure with density above saturation density will eventually stop, due to the transition to some other phase 
of matter, either previously considered scenarios such as pion and kaon condensation |5^, or superconducting quark 
matter. In neutron stars, K~ condensation is additionally favored because of a large electron chemical potential. 

Although lacking explicit confinement, the instanton model does provide the interaction to bind quarks into nucle- 
ons . The question whether the model provides the repulsive core necessary to prevent nucleons from collapsing into 
6 and more quark clusters remains open. Nevertheless, it seems plausible that a sufficiently sophisticated treatment 
could lead to nuclear matter as the correct ground state at small density. 

Can there be other phases, in addition to nuclear matter and superfluid quark matter? In the remainder of this 
section we will discuss a number of possibilities connected with diquark fluids or Bose condensates. These diquarks 



would not be Cooper pairs, but tightly bound states. As discussed in sect. [II, the instanton model seems to predict 
such states as bound scalar ud diquarks. The corresponding energy per baryon E/B = 3Mdq/2 ~ 800 — 900 MeV is 
lower than the nucleon mass, so it seems natural to look for a diquark phase. 

Naively, one would expect diquarks at T = to be Bose-condensed in the zero momentum state, because in 
addition to the gain in binding energy, there is a gain over nuclear or quark matter because no Fermi motion is 
required. However, since diquarks are colored, this phase could not be color neutral. This means that we have to 
consider either (i) significant motion of diquarks or (ii) add color- neutralizing quarks. 

Let us start with the first idea, with only diquarks at T = 0. A good starting point is the question why - if the scalar 
diquark is bound ~ a two-baryon state, such as the deutcron, does not decay into a more tightly bound three-diquark 
state. The simplest color singlet combination of three {ud)i — (j)i scalar diquark fields is 4'ii4'i24'a^ii,i2,i'ii but this 
wave function is antisymmetric, violating Bose statistics. This means that we have to consider p-wave diquark states 
(/)™ {e.g., in a bag), where m is the third component of angular momentum, in a symmetric combination: 

V = <Ptl <Pi2 fee eml,m2,m3- (111) 

Simple estimates show that such a state is no longer more economical. 
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Similarly, one can construct a wave function for infinite diquark matter starting from a set of plane wave states. 
The ground state would then be a quite peculiar "color crystal" , where the diquark momentum is determined from 
the balance of kinetic and potential energy, the latter resulting from the color-electric fields. We have not attempted 
to calculate the energy in this phase, but it seems much more natural to consider a quark-diquark phase. 

The quark-diquark (QDQ) phase could occur as an intermediate phase between nuclear matter and the color 
superconducting phase. Let us focus on a flavor composition corresponding to neutron matter (relevant for dense 
stars): ud diquarks plus an equal amount of d quarks. In this case, the total color and electric charge are zero. 
However, if the density is very low, color has to be neutralized over large distances, and confinement prohibits a phase 
like this. It is amusing to note that even without confinement, there is no low density QDQ phase. One reason is 
that the nucleon is bound w.r.t. a diquark and a quark. Another reason is that at very low density Fermi motion 
in the QDQ phase is more costly than in the nuclear phase. For example, let us set the threshold for the QDQ 
phase, Mdqq = M^q + M equal to the nucleon mass niN- Then the density of color-compensating d quarks is equal 
to the density of neutrons in the nuclear phase at the same baryon density. Furthermore, since both have the same 
degeneracy factor (gs=2 due to spin, the color being fixed by the color of the diquark condensate), both the d quark 
and the neutron have the same Fermi momentum pf. The kinetic energy pp/{2M) is then smaller for the neutron 
because of its larger mass. 

Nevertheless, it is not obvious which phase is preferred at densities a few times nuclear matter density. The QDQ 
phase is very different from both nuclear matter and color superconducting {Nf = 2) quark matter, in particular both 
chiral symmetry and color are broken. This suggests that if such a phase exits, it is probably an isolated minimum, 
separated by first order transitions on both sides. Let us try to estimate if such a window may exist. We assume 
that {qq) is large in this phase, still providing an effective quark mass 0(400 MeV), and therefore this phase should 
have no strange component, in contrast to the CSC3. The interactions in nuclear matter crucially combine long range 
attraction with short range repulsion. A similar description may be approximately valid for diquarks as well. Let us 
use a simple model, with only repulsive interactions represented by the scattering length|^ a. The energy per baryon 
of the diquark Bose gas is 



where Udq denotes the diquark density. The first term is just the mean-field interaction of the condensed diquarks, the 
second term stems from non-condensed bosons, as follows from the classic Lee- Yang paper js^ ]. In this approximation 
the unphysical behavior of the ideal Bose gas is overcome, the chemical potential grows with the density, and at 
some density it becomes favorable to split some diquarks into quarks. We are then led to a mixture of a Bose gas 
of diquarks and a Fermi gas of quarks (now of the same color), with chemical potentials related by the equilibrium 
condition = 2/x. Such a description leads to a more natural transition to quark matter with Cooper pairs at high 
density. 



^The radius of the nucleon repulsive core is about 0.4 fm, but diquarks (and constituent quarks) are smaller objects. If they 
are instanton-generated, their core should be of the order of the typical instanton radius p ~ 1/3 fm |^^, which we took as a 
representative value. 
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FIG. 14. Pressure versus baryonic chemical potential (3 times /i for quarks as used above) for 3 phases: nuclear matter, the 
quark-diquark phase made of a Bose gas of interacting diquarks and a Fermi gas of constituent quarks (no bag constant), and 
a Fermi gas of u, d, s quarks with current masses (for two different bag constants). 

So far we have ignored the role of confinement. One possibility is that the transition to the quark-diquark phase 
leads to deconfinement. If not, one has to include the energy associated with separated color charges. It is well known 
that the confining potential in vacuum is linear V{r) = Kr, with a string tension if ~ 1 GeV/fm. However, for small 
r this relation only holds for heavy, point-like quarks, not for light ones. Various non-relativistic models of hadronic 
structure use some effective K , reduced by a significant factor, in order to obtain a good description of hadronic 
masses. Indeed, it was was found on the lattice - by measuring the string tension after smoothing the gauge fieldfl - 
that the effective potential between constituent quarks is very small at small r, but approaches V{r) — Kr at large 
distances. The motivation for smoothening is the due to the extended nature of constituent quarks, as opposed to 
essentially point-like heavy flavors (c, 6). To put it differently: color-strings only form if constituent quarks do not 
overlap. For such a potential the average energy of colored strings seems to be negligible for the relevant densities 
considered above. 



C. Observable signatures 

In this section we would like to discuss potential experimental signatures of the quark-diquark and quark super- 
conducting phases, in particular with regard to heavy-ion experiments and neutron stars. Let us start with the 
quark-diquark phase. It is easy to estimate the critical temperature for the quark-diquark mixture by applying 
Einstein's ideal gas expression for Bose condensation 

T, = iMnJ^/Mdq . (113) 

Assuming ndq = 3*no and Mdq = 0.6 GeV we find Tc — 120 MeV, which is expected to be further reduced by the short 
range repulsive core. Heavy-ion collisions at SIS/BEVALAC energies (1-2 AGeV), where comparable compression is 
reached, lead to a heating of the system of up to T « 100 MeV, and so we conclude that even if this phase exists and 
our estimates are valid, it can be involved only peripherally. The same is even more true for the color superconducting 
phase. Because of the high density and low critical temperature of the color superconductor, it is not likely to be 
produced in heavy-ion collisions. 



^For example, the results of can be approximated as V = K\R — Ro\9{R — Ro) with the standard string tension K ~ 
1 GeV/fm but a rather large Ro — 0.7 fm. 
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Even if the gap is still larger than expected so that matter in the color superconducting phase would be produced in 
heavy-ion collisions, its presence will be hard to establish. The two most spectacular manifestations of superconduc- 
tivity, perfect conductivity and the Meissner effect, are very difficult to detect for a short-lived sample. The transition 
to the superconducting state has an effect on the equation of state, but the condensation energy e ^ //A^pi^/(27r^) 
is small compared to the energy density of a Fermi gas. In this context, the effect of quark-diquark phase is likely 
to be larger. Finally, the superconducting phase, in particular color-flavor-locking, may have an effect on the flavor 
composition. But this effect is not very specific since strangeness enhancement is a general consequence of quark 
matter formation. 

For these reasons, compact stars and stellar explosions are probably a more appropriate place to search for ob- 
servational consequences of quark superconductivity. For a recent review on the structure of neutron stars we refer 
the reader to Refs. [ ^J59| . Owing to theoretical uncertainties in the nuclear equation of state, the central density of 
neutron stars is not very well known. The main experimental constraint comes from the fact that neutron stars with 
masses M^s = 1-45M5„„ have definitely been observed. This still allows for central densities as low as 3no or as high 
as lOno- Ultimately, a better handle on the central density will come from measurements of neutron star radii. 

In neutron star structure calculations quark matter is almost always treated as a simple Fermi gas, confined by a 
bag constant B. Let us only note here that there are two distinct scenarios: (i) For sufficiently small values of B 
strange quark matter is absolutely stable, and the entire star is in this phase, while (ii) for large B the outer part 
of the star consists of nuclear matter, while the interior contains various mixed phases (quark matter sheets, rods or 
clusters). These funny arrangements owe their existence to the possibility to move charge from the quark phase to 
the hadronic phase. The shape is then determined by the interplay of the equation of state and long range Coulomb 
forces. Quark superconductivity will again have an influence on the equation of state, but as before we expect this to 
be a correction on the order of 0((A/^)^), small compared to the uncertainties in the bag constant. 

In neutron stars, a more direct measure of the gap is provided by the cooling history of the star. Without a 
gap, neutron stars can efficiently cool by /3-decay of thermally excited u, d quarks. This process seems to lead to 
unacceptably large cooling rates Jsot . Such a problem emerged already for neutrino emission from nuclear matter. 
In that case it is solved due to nuclear superfluidity. Since both neutrons and protons are gapped, there are no 
single-particle states in the vicinity of the Fermi surface. For two-flavor quark matter, a similar problem arises since 
the quarks of the third color have no or only very small gaps. It is absent in the color-flavor-locked phase, because 
all quarks acquire a gap. For realistic values of the strange quark mass, the system is likely to be in the two-flavor 
phase at moderate densities, so that the "cooling-problem" for neutron stars might still persist. 

Another feature of the CSC-phases is that Cooper pairs are electrically charged, so one might expect the color 
superconductor to be electrically superconducting as well. This is not quite true. For both the Nf = 2 and color- 
flavor locked phase, there is a modified charge operator that is not broken, so there is a linear combination of the 
photon and the diagonal gluons that remains massless. Nevertheless, since the photon inside the superconducting 
phase is different from the photon outside, magnetic flux will be partially expelled. Magnetic fields in pulsars are very 
large, up to 10^^ Gauss. Fields on the order of 10^^ Gauss were recently suggested to drive the so called "magnetars" , 
and 10^^ gauss is the absolute upper limit allowed by star stability (virial theorem). However, even such fields are not 
yet large enough to significantly influence the CSC phase. On the other hand, simple scaling considerations show that 
in order to have the field completely expelled from the quark matter core would cost energy of order 0{R^) (where R 
is the star radius), while transferring the field lines through some channels into the superconductor only costs 0{R). 
Inside these channels the field should be at the critical value Be ^ 10^^ Gauss. They can be either macroscopic (if 
the superconductor is of the first kind) or microscopic Abrikosov vortices (if it is of the second kind) . 

X. SUMMARY AND CONCLUSIONS 

We have studied the interplay of instantons, superfiuidity /-conductivity and chiral symmetry breaking in QCD at 
finite density. Unlike many schematic models based on short-range interactions abstracted from one gluon exchange, 
the instanton model has the virtue of providing a realistic phenomenology of the zero-density ground state, including 
such features as spontaneous chiral symmetry breaking and a correct description of (light) hadron spectroscopy. This 
gives us some confidence for a semi-quantitative investigation of cold quark matter at moderate densities as might be 
encountered, e.g., in the core of neutron stars. 

We started out by reviewing some properties of the instanton-induced quark-quark interaction at /i = 0: it originates 
from the same effective interaction that leads to chiral condensation and a light pion in the QCD vacuum, predomi- 
nantly acting in the scalar-isoscalar, color antitriplet qq channel. We found these correlations to be sufficiently strong 
to generate a bound-state pole in the corresponding diquark propagator. 

In a second step we studied the density dependence of the instanton-induced ('t Hooft) interaction. It arises through 
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the modification of the quark zero modes as we introduce a chemical potential. Both the density-dependence of the 
instanton form factors, governing the effective quark- (anti-) quark interactions, and of the instanton-antiinstanton 
overlap matrix elements, relevant for a statistical treatment of the instanton liquid, have been assessed. The main 
difference to earlier calculations (using local approximations or schematic formfactors) manifests itself in quantitatively 
somewhat larger superconducting gaps (~ 150-200 MeV). 

In the first main part of this article we have performed a systematic study of the finite-/^, = 3-QCD phase 
structure for 2, 3 and 2-1-1 flavors within the mean- field approximation: 

(i) For two massless flavors we confirmed the usual transition to a ud superconductor at /Xc — 300 MeV with gaps 
A ~ 200 MeV. Apart from some quantitative deviations the use of microscopic instanton formfactors entailed 
the possible existence of a novel intermediate phase, namely chirally broken quark-diquark matter (characterized 
by simultaneous chiral symmetry breaking and diquark condensation). 

(ii) For three massless flavors we returned to a simplified treatment based on a idealized sharp formfactor (amenable 
to the more involved calculations with somewhat less accuracy). We found that, like schematic one-gluon 
exchange interactions, instantons lead to color-flavor locking, despite their very different color-flavor vertex 
structure. In addition, there emerged an interesting new feature: in contrast to OGE, instantons generate a 
non-zero (albeit small) chiral condensate in the superfluid phase. 

(iii) The consequences of finite (current) strange quark masses rus have been investigated. With increasing m,s a 
sharp phase transition from the color-flavor locked phase to the two-flavor siiperconductor occurs. At small 
chemical potential, the critical strange quark mass is smaller than the physical mass. As the chemical potential 
grows, the critical mass is also expected to grow. We find that the color-flavor locked phase will appear only 
at /i > 450 — 500 MeV. This implies that chiral symmetry would be restored initially and then broken again at 
higher density, but with much smaller condensates. 

In the second main part of our article a statistical mechanics treatment of the instanton liquid has been employed. 
This enabled us to incorporate effects that go beyond the standard mean-field approximation, in particular those 
associated with instanton-antiinstanton molecules. We demonstrated how to handle a complex fermion determinant 
in this context (using a gluonic interaction that has been averaged over the relative color orientation). Without 
molecules, the results were shown to be in reasonable agreement with the mean- field analysis in the first part. Including 
correlations, molecule formation constitutes a ^10% effect in the free energy /critical chemical potential of the T=0-, 
Nf = 2-transition. However, with rising temperature I-A molecules are expected to play an increasingly important 
role in the transition between the superconducting and the plasma phases, corresponding to the c;ritical temperature 
for superconductivity. Furthermore, the statistical mechanics approach allows to assess the /i-dependence of the total 
instanton density. It turns out that the latter is indeed approximately constant for all /U-values under consideration, 
which is not really surprising as the gluonic energy within the instanton component carries the dominant fraction of 
the free energy of the system. This supports the respective assumption made in earlier works as well as in the first 
part of this article. 

Based on our results, we conjectured a qualitative picture of the QCD phase diagram in the T-/i plane. We 
argued that for realistic values of the strange quark mass, there are both finite-/Lt and finite-T transitions between 
the color-flavor locked (CSC3) and two-flavor superconducting (CSC2) states. With increasing density and at small 
T, chiral symmetry is first restored and then (weakly) broken again. We also elucidated on the possibility of a new 
quark-diquark phase, in which nucleons are dissolved into a Bose gas of wd-diquarks and a Fermi gas of unpaired 
quarks, characterized by simultaneously broken color and chiral symmetry. It represents a natural possibility for an 
intermediate phase between the hadronic and superconducting phase, although numerical estimates for its existence 
are rather uncertain. 

Concerning experimental c;onsequcnc;cs. we conclude that heavy-ion reactions are unlikely to reach into the rich 
high-density/low- temperature phase structure discussed here. In this respect neutron stars are much more promising. 
The pairing gaps should leave their traces in cooling rates and even in the equation of state provided the pairing 
gaps are large enough. Decomposition of the "external" into "internal" magnetic fields presumably imply complicated 
configurations inside the stars. 

As an outlook, we expect that it should be rather straightforward to generalize our approach to a simultaneous 
account of finite /i and T. Another major challenge is to go beyond the mean-field approximation in the quark sector 
in order to address clustering of quarks into nucleons, and, even more difficult, the nucleon-nucleon interaction. The 
instanton liquid model does provide the correct properties on the one-nucleon level. Turning the instanton liquid 
into a realistic description of nuclear matter, however, requires a long way to go. These and related issues certainly 
provide exciting opportunities for future research in the field of finite-density QCD. 
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APPENDIX A: FIERZ TRANSFORMATIONS 



Let us denote a general four-fermion interaction of the type (^) by {rpaOipb){'4'cOipd) ■ An effective "mesonic" 
interaction corresponds to the sum of the direct and exchange {s+u) channels. The corresponding kernel is the sum 
of the original interaction and its Fierz transform: 



(Al) 



An effective diquark (f-channcl) interaction is obtained by first transposing and then performing the Fierz transfor- 
mation 



(A2) 



We use Euclidean gamma matrices in the standard or chiral representations, C = — 172707 C*^ — C^C — 1. If the 
isospin and color parts of O are 1, e = — 1 if the Dirac part of O is 7^, (t^,^, and e = 1 when it is l^, 75, j^j^. For the 
flavor part, rj = —T2 and t^^ = ri,3. 

For the Fierz transformation the following completeness relations are used: 



for color, 



5aa'5bb' 



= -^Sab'S 
= t^Sab'S 



BA' 



'''AB''''BA' 



BA' 



:'''ab'^ba' 
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(A4) 
(A5) 



2 2 

for isospin, and, defining S" = 1 (g) 11,^ = 7^ (g) 7^,, P = 75 ® 75, A = 757^ ® j^,T = cr^i, (g) cr^i,, with cr^^ = 1/2[7m' 7i^]: 



/S\' /1/4 1/4 -1/8 -1/4 l/4\ 
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(A6) 



for the Dirac structures. Due to the tensor product structure of the matrices O, one has to Fierz transform separately 
the Dirac, color and isospin parts by using the above relations and multiply them accounting for the fact that fermion 
fields anticommute. The latter generates an additional (-1) in the w-channel, and for the t-channel it excludes any 
symmetric parts of the O' matrices. The final results are given in the main text, Eqs. (Q0[) and (pT|). 



APPENDIX B: INSTANTON FORM FACTORS AT FINITE ^ 

The two structures in the form factor are 

1 r/ e 

1 - 2— — - - fit 



B = Api / dRR^ / dri sir? n ^= 

Jo Jo R^/W 



R^+p^ 



sm(pr) 



- 2 



t 



R^ + p'~ 



p I cos{pr) 



sm{kr) 
kr 



(Bl) 



-4p / dRR'^ / di] sin ry 
/o Jo 
' pt 2r 



R./WT7 



2t 



cos{pr) 



cos(fcr) sin(fcr) 
kr k'^r^ 



R^+p^ 

^ — ILUt 



sin (/ir) 



(B2) 



Let us separate the real and imaginary parts of A and B: A ^ Ac + iAs, B ~ Bs + iBc, the subscript referring to 
cos{ut) or sin{ujt). We also introduce k^ = k ± p. Then 
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k 



dR / dr] 



2R sin 7] 



R sin Tj cos{ujR cos rj) 



D2 , 2 + 7^^^ — ) cos(A:+i?sin7?) 
R-' + p-' kR smt] J 

+ [ 1^- ,^/r>9^. „9^ ) sm{k+R sin r]) + (/i ^ -/x) 



/ dii / 

Jo Jo 



dr] 



R cos T] sin(wi? cos r]) 



1 2 \ 

TT:^ 9 1" 771^5 ^ ^ M cos(fc+i?sin?7) 

kR^sin^r] k{R^ + p^) J 

fc+ 2i?sin?7\ . , „ . , 
+ ( T-PT^ 1- t;;^— — ^ I sin(A;+i?sin7y) 



Br, 



kR sin 7/ R? + p"^ 

cos{ujRcosr]) 



dr]- 



2i?2 cos^ 



cos(fc"'"iisinj7) + /ui?sinr/sin(fc"'"i?sinr/)] — (/i — > — /x) 



dR \ dr] 



R cos T] sin(a;i? cos r]) 



, „ . , 2i?sin77 N 
cos(«;"^ it sin ?7) — — 5- sin(fc"^ it sm 77) 



i?2 + p2 



Using the following basic integrals, 



d?? cos (a cos ry) cos (/3 sin 77) = 7rJo(\/ o? + /J^) 



and, denoting = -^/w^ -)- (fc='=)2, 
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one can rewrite the above expressions as derivatives of I^K^ according to 
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' integral over /i' is just ^io^ifj^. Furthermore, 
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Denoting £4K+ = I+K+ - 1+ K+ = A+, we get 
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These Fourier transforms of the fermion zero modes at finite /i agree with the results obtained in [p^. 

APPENDIX C: GRAND CANONICAL POTENTIAL IN THE CORNWALL-JACKIW-TOMBOULIS 

(CJT) FORMALISM 

To derive the grand canonical potential in MFA we start from a generating functional with bilocal meson and 
diquark sources, 

eiip{W[J,J,K,K])= [ V^V^P exp \ (So + U + [ ipJ^p + [ Jip^ + [ iPKi;'^ + [ i;'^K^P)\ , (CI) 



where So is the free fermion action and U = J {£mes + ^diq) with £mes,£diq given in Eqs. ( |lO| , [ri|) . Using a matrix 
representation, we write the free part of the action in momentum space as 



5o 



where the sources J, J, K have the following properties: 

- r{-p) = j{p) 

-K^i-p) = K{p) 
~K^{-p) = K{p) . (C3) 



The path integral over the fermion fields is 



e'^" = / I?V2?^e^ = det5 ( _ f ^ ^^^o ) = det 



T4^o = itrlnA^ , (C4) 

where the trace includes the momentum integration, and the products between the fields (sources) include convolutions 
in momentum as well as all other indices. Using the identities 

B a) ^ yB ) yB-^A 1 
Indet I = Indet(-BB) + trln(l + (^^^i^ 

= lndet(-i?B) + tr)^— — 

n— 1 ^ 



Indet(-BB) + tr ^ - — (B-^AB-^A + B^^AB^^A) 



2n 

n=l 

= Indet(-BB) + trln(l - B^^AB^^A) 

= tr ln(-BB + BAB^^ A) , (C5) 

Wo - ^trln[-(J + i)(J - ^00^^) + (J + ^Gq i)K(J + IGoY'^ ■ (C6) 

Next, one introduces the classical two-point fields F = F = G = G = ^j- and performs a Legendre 
transformation 



we obtain 
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r[G, G, F, F] = W[J, J, K, K] - tr(GJ + GJ + KF + KF) (C7) 
with J, J, K, K expressed as functionals of G, G, F, F. From the non-interacting part Wq we have 

G - ^( J - 1^ - X( J + Igo 
F = -(J+iGo-i)-iifG 

F = -(J-iGo-i^)-iif)-ii^G, (C8) 

and 

Fo = -itrln(-GG + GFG-^F) + ^tr(Go ^G - Gq ^"^G - 2) . (C9) 
Following the CJT approach ||3^] for the four-fermion interaction, one can show that 

where Sint is the original four-fermion interaction and we have explicitly indicated the dependence on Ti. Now it 
becomes clear that, in order to perform the functional derivatives, it is conveni ent t o Fierz-rearrange Sint into the 3 



channelss, t and u as outlined above. The sum of the four derivatives in Eq. (CIC) naturally suggests the use of U 
in Eq. ( |Clf ) with two- by- two derivatives w.r.t. the fermion fields that are contracted with the same matrix O*. 

The Hartree-Fock scheme is equivalent to (i) neglecting F4 , which is the sum of all four-particle irreducible diagrams 
and is of order 0{Ti^), and (ii) considering only translationally invariant solutions {e.g., G{p,p') = G{p)S^{p — p'), 
and the same for F, F). The mesonic and the diquark terms in the above lowest-order interaction term have the form 
1/4(M,)(M,), 1^4(D,)(D,), where 

D. = j ^tr(Ftgf'^.FTo'GF*) , (Gil) 

i.e., they are products of mesonic and diquark condensates. A priori one might consider condensation in any channel, 
and then study if it is energetically favored. It is reasonable to start with the channels that are most attractive 
[e.g., those which generate bound states in vacuum). As we have seen in sects. |l| and these are the scalar color 
singlet mesonic channel and the diquark color-3 one. For the Nf — 2 case there is only one pattern of diquark 
condensation, which necessarily breaks the color symmetry down from S'[/(3)c — > SU(2)c. The diquark condensate 
is proportional to a unit vector in SU{i)c. Moreover, because the flavor part of the vertex is antisymmetric (due to 
T2), the Pauli principle requires that the above unit vector belongs to the antisymmetric part of SU{'i)c. Without 
loss of generality we can choose this vector to be the Gell-Man matrix A2 (recall that all Gell-Mann matrices in this 
paper are normalized to 3, i.e., trA? = 3). We can identify the unbroken SU{2)c with the upper 2x2 corner of the 
SU{i)c group. If both qq and qq condensates are present, the latter should consist of two parts - one that involves 
the two colors from the unbroken SU{2)c and the second one which involves quarks and antiquarks from the third 
color. One should note that there are two terms in Cmes contributing to the qq condensates: the isoscalar color singlet 
one and the isoscalar SU{3)c octet one proportional to Ag. Both have projections onto the SU{2)c singlet as well 
as onto the subgroup represented by third color of (unpaired) quarks. Owing to the two different chiral condensates, 
it is convenient to split the propagators G, G into a SU{2)c part Gi,Gi and a part involving the third color only, 
G2 , G2 . The F and F propagators are proportional to the diquark condensate and hence to A2 . Putting everything 
together, Eq. ( C1C| ) becomes 



r = -itrln(-GiGi + GiFGf^F) -f itr(Go ^Gi - ^^Gi - 2) 
-itrln(-G2G) + ^tr(G^iG2 - G^'^^G2 - 2) 



41 



+5g^^F2(]J-:Y)tr(FC75A2r2/3)tr(/3*C75A2T2F), (C12) 



where a and /? are the formfactors defined in Eqs. (p7D,(|3£ 

Since G, G, F, F are classical fields, F should vanish under the respective variations. 

APPENDIX D: GAP-INDUCED DAMPING OF SINGLE-QUARK PROPAGATORS 

In the superconducting phases, single-quark propagation in the vicinity of the Fermi surface is damped in the 
temporal direction due to the presence of the finite energy gap. This effect has to be included in the evaluation of 
the I- A overlap matrix elements, Eq. (^o|), at finite density (this is particularly crucial for the activity Zd of Eq. (|9l|), 
which otherwise would diverge). Rewriting the overlap matrix elements as 

TjA = - I d'^x [4>\{x - zi\ ^ - 1^74)] {i irn^y^ [(« ?> - «Ai74)*o,A(a; - za; ^JL)\ , (Dl) 



they are readily interpreted as a quark hopping amplitude, represented by two amputated (anti-) instanton vertices 
and an intermediate (free) quark propagator (i ^ — ^^74)"^. In the following we will evaluate an approximate damping 
factor for this propagator. Starting from the expression of the (massless) momentum space propagator at finite /i in 
the superconducting phase (ignoring the Dirac structure) 



G(p,M,A) = -2 — ^ — -2 



Po ~ — + irj 
ul 



(D2) 



po - e(p) Po + e(p) - ir] 

— Up — fi, e{pY = + A^, Mp = i[l + ^/e(p)], = ^[1 — ^/e(p)]), we compute its Fourier transform as 

G(z;a.,A)= j ^e''P^G{p,t,,A) 

00 

L_ f pdp sin(pr) e-'^(f)^* . (D3) 



For A ^ one recovers the result | 
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G{z;fj,) 
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(2z4 + ^z^) cos{^ir) + (zl — + ^z^z^) ^^^^^^^^ 

r 



27r2z4 

An approximate correction to T7^(z; /i) is thus obtained by supplying it with the ratio 



(D4) 



R{z4; n, A) = — — — , (D5) 

g(z4;m) 

where we have restricted the l.h.s of Eq. ( ^5| ) to r=0 to avoid artificial singularities caused by oscillations in G(z4, r; /i). 
However, since (for Nf = 2) only two out of three quarks at the Fermi surface can participate in the diquark condensate. 



the correction (D5) enters on average with a smaller power, i.e., 

Tja{z; A) ~ R{zr, fi, A)2/3 Tja{z; fi) . (D6) 

The net effect of R{zi] /i. A) is a damping of the zero mode propagation, which, in fact, is much less pronounced than 
the naive expectation, oc e"'^^*, would suggest. On the other hand, ilinst disfavors finite values for A. 
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